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Abstract 

We propose a review of recent developments on entanglement and 
non-classical effects in collective two-atom systems and present a uni- 
form physical picture of the many predicted phenomena. The col- 
lective effects have brought into sharp focus some of the most basic 
features of quantum theory, such as nonclassical states of light and en- 
tangled states of multiatom systems. The entangled states are linear 
superpositions of the internal states of the system which cannot be 
separated into product states of the individual atoms. This property 
is recognized as entirely quantum-mechanical effect and have played 
a crucial role in many discussions of the nature of quantum measure- 
ments and, in particular, in the developments of quantum commu- 
nications. Much of the fundamental interest in entangled states is 
connected with its practical application ranging from quantum com- 
putation, information processing, cryptography, and interferometry to 
atomic spectroscopy. 
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1 Introduction 



A central topic in the current studies of collective effects in multi-atom sys- 
tems are the theoretical investigations and experimental implementation of 
entangled states to quantum computation and quantum information pro- 
cessing jT]. The term entanglement, one of the most intriguing properties of 
multiparticle systems, was introduced by Schrodinger in his discussions of 
the foundations of quantum mechanics. It describes a multiparticle system 
which has the astonishing property that the results of a measurement on one 
particle cannot be specified independently of the results of measurements on 
the other particles. In recent years, entanglement has become of interest not 
only for the basic understanding of quantum mechanics, but also because 
it lies at the heart of many new applications ranging from quantum infor- 
mation Ej , cryptography |2] and quantum computation [HI E] to atomic 
and molecular spectroscopy jHUHl- These practical implementations all stem 
from the realization that we may control and manipulate quantum systems 
at the level of single atoms and photons to store and transfer information in 
a controlled way and with high fidelity. 

All the implementations of entangled atoms must contend with the con- 
flict inherent to open systems. Entangling operations on atoms must provide 
strong coherent coupling between the atoms, while shielding the atoms from 
the environment in order to make the effect of decoherence and dissipation 
negligible. The difficulty of isolating the atoms from the environment is the 
main obstacle inhibiting practical applications of entangled states. The envi- 
ronment consists of a continuum of electromagnetic field modes surrounding 
the atoms. This gives rise to decoherence that leads to the loss of informa- 
tion stored in the system. However, it has been recognised that the collective 
properties of multi-atom systems can alter spontaneous emission compared 
with the single atom case. As it was first pointed out by Dicke [lOj, the inter- 
action between the atomic dipoles could cause the multiatom system to decay 
with two significantly different, one enhanced and the other reduced, sponta- 
neous emission rates. The presence of the reduced spontaneous emission rate 
induces a reduction of the linewidth of the spectrum of spontaneous emis- 
sion fUE]. This reduced (subradiant) spontaneous emission implies that 
the multi-atom system can decohere slower compared with the decoherence 
of individual atoms. 

Several physical realisations of entangled atoms have been proposed in- 
volving trapping and cooling of a small number of ions or neutral atoms ^3 
ITi| US] • This is the case with the lifetime of the superradiant and sub- 
radiant states that have been demonstrated experimentally with two barium 
ions confined in a spherical Paul trap [T^ IT^ . The reason for using cold 
trapped atoms or ions is twofold. On the one hand, it has been realised that 
the trapped atoms are essentially motionless and lie at a known and con- 
trollable distance from one another, permitting qualitatively new studies of 
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interatomic interactions not accessible in a gas cell or an atomic beam [T7] . 
The advantage of the trapped atoms is that it allows to separate collective 
effects, arising from the correlations between the atoms, from the single-atom 
effects. On the other hand it was discovered that cold trapped atoms can be 
prepared in maximally entangled states that are isolated from its environ- 
ment dHl [111201 ini 1221 ■ 

An example of maximally entangled states in a two-atom system are 
the superradiant and subradiant states, which correspond to the symmetric 
and antisymmetric combinations of the atomic dipole moments, respectively. 
These states are created by the interaction between the atoms and are charac- 
terized by different spontaneous decay rates that the symmetric state decays 
with an enhanced, whereas the antisymmetric state decays with a reduced 
spontaneous emission rate. The reduced spontaneous emission rate of the 
antisymmetric state implies that the state is weakly coupled to the envi- 
ronment. For the case of the atoms confined into the region much smaller 
than the optical wavelength (Dicke model), the antisymmetric state is com- 
pletely decoupled from the environment, and therefore can be regarded as a 
decoherence-free state. 

Another particularly interesting entangled states of a two-atom system 
are two-photon entangled states that are superpositions of only those states of 
the two-atom system in which both or neither of the atoms is excited. These 
states have been known for a long time as pairwise atomic states or multi- 
atom squeezed states [23 1211 123 123 123 12H1- The two-photon entangled states 
cannot be generated by a coherent laser field coupled to the atomic dipole 
moments. The states can be created by a two-photon excitation process 
with nonclassical correlations that can transfer the population from the two- 
atom ground state to the upper state without populating the intermediate 
one-photon states. An obvious candidate for the creation of the two-photon 
entangled states is a broadband squeezed vacuum field which is characterised 
by strong nonclassical two-photon correlations [^ I3UI [M] . 

One of the fundamental interests in collective atomic effects is to demon- 
strate creation of entanglement on systems containing only two atoms. A 
significant body of work on preparation of a two-atom system in an entan- 
gled state has accumulated, and two-atom entangled states have already been 
demonstrated experimentally using ultra cold trapped ions in free space |14| 
122] and cavity quantum electrodynamics (QED) schemes [S3 [Ml- In the free 
space situation, the collective effects arise from the interaction between the 
atoms through the vacuum field that the electromagnetic field produced by 
one of the atoms influences the dipole moment of the another atom. This 
leads to an additional damping and a shift of the atomic levels that both 
depend on the interatomic separation. In the cavity QED scheme, the atoms 
interact through the cavity mode and in a good cavity limit, photons emit- 
ted by one of the atoms are almost immediately absorbed by the another 
atom. In this case, the system behaves like the Dicke model. Moreover, the 
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strong coupling of the atoms to the cavity mode prevents the atoms to emit 
photons to the vacuum modes different from the cavity mode that reduces 
decoherence. 

Recently, the preparation of correlated superposition states in multi-atom 
system has been performed using a quantum nondemolition (QND) measure- 
ment technique |35j. Osnaghi et al. [36j have demonstrated coherent control 
of two Rydberg atoms in a non-resonant cavity environment. By adjust- 
ing the atom-cavity detuning, the final entangled state could be controlled, 
opening the door to complex entanglement manipulations [SZl- Several pro- 
posals have also been made for entangling atoms trapped in distant cavi- 
ties |SH1 EHl UHl El 113 113, or in a Bose-Einstein condensate jH |1S|. In a 
very important experiment, Schlosser et al. jlHI succeeded in confining single 
atoms in microscopic traps, thus enhancing the possibility of further progress 
in entanglement and quantum engineering. 

This review is concerned primarily with two-atom systems, since it is gen- 
erally believed that entanglement of only two microscopic quantum systems 
(two qubits) is essential to implement quantum protocols such as quantum 
computation. Some description of the theoretical tools required for predic- 
tion of entanglement in atomic systems is appropriate. Thus, we propose to 
begin the review with an overview of the mathematical apparatus necessary 
for describing the interaction of atoms with the electromagnetic field. We 
will present the master equation technique and, in addition, we also describe 
a more general formalism based on the quantum jump approach. We review 
theoretical and experimental schemes proposed for the preparation of two 
two-level atoms in an entangled state. We will also relate the atomic entan- 
glement to nonclassical effects such as photon antibunching, squeezing and 
sub-Poissonian photon statistics. In particular, we consider different schemes 
of generation of entangled and nonclassical states of two identical as well as 
nonidentical atoms. The cases of maximally and non-maximally entangled 
states will be considered and methods of detecting of particular entangled 
and nonclassical state of two-atom systems are discussed. Next, we will ex- 
amine methods of preparation of a two-atom system in two-photon entangled 
states. Finally, we will discuss methods of mapping of the entanglement of 
light on atoms involving collective atomic interactions and squeezing of the 
atomic dipole fluctuations. 

2 Time evolution of a collective atomic sys- 
tem 

The standard formalism for the calculations of the time evolution and correla- 
tion properties of a collective system of atoms is the master equation method. 
In this approach, the dynamics are studied in terms of the reduced density 
operator pA of the atomic system interacting with the quantized electromag- 
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netic (EM) field regarded as a reservoir j^TjEHlllSl- There are many possible 
realizations of reservoirs. The typical reservoir to which atomic systems are 
coupled is the quantized three-dimensional multimode field. The reservoir 
can be modelled as a vacuum field whose the modes are in ordinary vacuum 
states, or in thermal states, or even in squeezed vacuum states. The major 
advantage of the master equation is that it allows us to consider the evolution 
of the atoms plus field system entirely in terms of average values of atomic 
operators. We can derive equations of motion for expectation values of an 
arbitrary combination of the atomic operators, and solve these equations for 
time-dependent averages or the steady-state. Another method is the quan- 
tum jump approach. This is based on the theory of quantum trajectories [20], 
which is equivalent to the Monte Carlo wave-function approach |5H , and 
allows to predict all possible trajectories of a single quantum system which 
stochastically emits photons. Both methods, the master equation and quan- 
tum jumps approaches lead to the same final results of the dynamics of an 
atomic system, and are widely used in quantum optics. 



2.1 Master equation approach 

We first give an outline of the derivation of the master equation of a system 
of N non-identical nonoverlapping atoms coupled to the quantized three- 
dimensional EM field. This derivation is a generalisation of the master equa- 
tion technique, introduced by Lehmberg jlTj, to the case of non-identical 
atoms interacting with a squeezed vacuum field. Useful references on the 
derivation of the master equation of an atomic system coupled to an ordi- 
nary vacuum are the books of Louisell j48j and Agarwal ||49j. The atoms 
are modelled as two-level systems, with excited state \ei), ground state \gi), 
transition frequency Ui, and transition dipole moments /ij. We assume that 
the atoms are located at different points ri, . . have different transition 
frequencies ui ^ 002 7^ ... 7^ uj^, and different transition dipole moments 

1^1 7^ f^2 7^ ■ ■ ■ ^ f^N- 

In the electric dipole approximation, the total Hamiltonian of the com- 
bined system, the atoms plus the EM field, is given by 

N 



/ 1 



«=1 ks 
N 

- ^^J2J2[f^i-9ksi^^{St + S~)a^^-ii.c.] , (1) 

ks «=1 

where = \ei){gi\ and S~ = \gi){ei\ are the dipole raising and lowering 
operators, 5*^ = (|ej)(ej| — \gi){gi\) /2 is the energy operator of the ith atom. 



which has wave vector k, frequency cj^ and the index of polarization s. The 



and dt^ are the annihilation and creation operators of the field mode ks, 
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coupling constant 



9ks I'"* 



ik-fi 



(2) 



is the mode function of the three-dimensional vacuum field, evaluated at the 
position Ti of the ith atom, V is the normalization volume, and e^^ is the 
unit polarization vector of the field. 

The atomic dipole operators, appearing in Eq. (P), satisfy the well-known 
commutation and anticommutation relations 
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(3) 



with (^f ) = 0. 

While this is straightforward, it is often the case that it is simpler to 
work in the interaction picture in which the Hamiltonian evolves in time 
according to the interaction with the vacuum field. Therefore, we write the 
total Hamiltonian (P) as 



H 



Ho + Hj 



where 



»=i ks ^ 

is the Hamiltonian of the non-interacting atoms and the EM field, and 

N 

Hi = -i^Y^Yl [f^i ■ 9ks (^0 {S^ + S^) hs - H-c. 



(4) 



(5) 



ks 



i=l 



(6) 



is the interaction Hamiltonian between the atoms and the EM field. 
The Hamiltonian Hq transforms the total Hamiltonian into 



H{t) 



JHot/h 



(h - Ho) e"'^"*/^ = V (t) 



(7) 



where 



V{t) = -^nJ:T.{f^^■9ksir^)Sta^ks^-'^-''--'^' 

ks «=1 

+/^.-^A;Jrl)5-a,.,e-^(-'^+-)*-H.c.} . 



(8) 



We will consider the time evolution of the collection of atoms interacting 
with the vacuum field in terms of the density operator paf characterizing the 
statistical state of the combined system of the atoms and the vacuum field. 
The time evolution of the density operator of the combined system obeys the 
equation 

a „ 1 r- 



dt 



Paf 



ih 



H, Paf 



(9) 



6 



Transforming Eq. Q into the interaction picture with 



p^^(t)=e^^°*/VAFe-^^°*/^ (10) 
we find that the transformed density operator satisfies the equation 



g-^AF it) 



1 

ih 



(11) 



where the interaction Hamiltonian V (t) is given in Eq. (jH)). 

Equation (jllj) is a simple differential equation which can be solved by the 
iteration method. For the initial time t = 0, the integration of Eq. leads 
to the following first-order solution in V (t) : 



Paf (t) = Paf (0) + ^ /„ dt' V (f) , p^p it') 



(12) 



Substituting Eq. (fT^ into the right side of Eq. (fTTj) and taking the trace 
over the vacuum field variables, we find that to the second order in V (t) the 
reduced density operator of the atomic system pA (t) = Tip Paf (^) satisfies 
the integro-differential equation 



^PA(t) = ^TTF[V{t),~PAFiO) 



1 



dt'TTF{[V{t),[V{t'),~pAFit')]]} . 



(13) 



We choose an initial state with no correlations between the atomic system 
and the vacuum field, which allows us to factorize the initial density operator 
of the combined system as 



~Paf (0) = Pa (O)pf(O) 



(14) 



where is the density operator of the vacuum field. 

We now employ the Born approximation j3H]; in which the interaction 
between the atomic system and the field is supposed to be weak, and there is 
no the back reaction effect of the atoms on the field. In this approximation 
the state of the vacuum field does not change in time, and we can write the 
density operator p^p (t'), appearing in Eq. (fT^. as 



~pAFit')=PAit')pFiO) . 



(15) 



Under this approximation, and after changing the time variable to t' = t — r, 
Eq. (|T!^ simplifies to 

^p(t) = ^Trp [f (t) , p (0) (Oy 

- 1 drTrp { [V (t) ,[v(t-T),p{t-T) ^ (0)]] } , (16) 
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where we use a shorter notation p = pA- 

Substituting the exphcit form of V {t) into Eq. (fTBj). we find that the 
evolution of the density operator depends on the first and second order cor- 
relation functions of the vacuum field operators. We assume that a part of 
the vacuum modes is in a squeezed vacuum state for which the correlation 
functions are given by [23 IHOl IHIj 



Tr 



Tr^ 
Tr^ 
Tr^ 
Tr^ 



pF (0) a^^ 

PF (0) Hfil^, 
PF (0) a\a^:,^, 
pF (0) flfc.a^,^, 



Tr? 



PF (0) 



ks 







5' {k 



k'] 5. 



k'] 5, 



D ituk)r N (tuk) + 1 
\DMfN{uJk)6^ {k 
{uJk)M{uJk)6' {2k. 
(iUk) M* (iUk) 6^ (2ks -k-k') 6s 



k-k']6 



(17) 



where the parameters (uk) and M (uk) characterize squeezing in the vac- 
uum field, such that N {uJk) is the number of photons in the mode k, M [uk) = 
\M (uk) \exp{i(j)s) is the magnitude of two-photon correlations between the 
vacuum modes, and (f)s is the phase of the squeezed field. The two-photon 
correlations are symmetric about the squeezing carrier frequency 2u!s, i.e. 
M (uk) = M {2ujs — UJk), and are related by the inequality 



\M < N (cok) (N i2us -cok) 



(18) 



where the term -fl on the right-hand side arises from the quantum nature of 
the squeezed field [201 EI] . Such a field is often called a quantum squeezed 
field. For a classical analogue of squeezed field the two-photon correlations 
are given by the inequality \M {ujk)\ < N (uk). Thus, two-photon correlations 
with < \M {uk)\ < N (uk) may be generated by a classical field, whereas 
correlations with (uk) < \M (cj^)! < (cuk) (N (2uJs — tOk) + 1) can only 
be generated by a quantum field which has no classical analog. 

The parameter D{ujk), appearing in Eq. ()17|) . determines the matching 
of the squeezed modes to the three-dimensional vacuum modes surrounding 
the atoms, and contains both the amplitude and phase coupling. The ex- 
plicit form of D (uk) depends on the method of propagation and focusing 
the squeezed field j2El ES]- For perfect matching, \D {ujk)^ = 1, whereas 
\D {uJk)f < 1 for an imperfect matching. The perfect matching is an ideal- 
ization as it is practically impossible to achieve perfect matching in present 
experiments |541 IH^ . In order to avoid the experimental difficulties, cavity 
situations have been suggested. In this case, the parameter D (uk) is identi- 
fied as the cavity transfer function, the absolute value square of which is the 
Airy function of the cavity j^HlEZI- The function \D (uk) P exhibits a sharp 
peak centred at the cavity axis and all the cavity modes are contained in a 
small solid angle around this central mode. By squeezing of these modes we 
can achieve perfect matching between the squeezed field and the atoms. In 
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a realistic experimental situation the input squeezed modes have a Gaussian 
profile for which the parameter D (uk) is given by [HZl US] 



D (uk) = exp —Wq sin^ 6^ — ikzj cos 6'^ 



(19) 



where 9k is an angle over which the squeezed mode k is propagated, and Wo 
is the beam spot size at the focal point Zf. Thus, even in the cavity situation, 
perfect matching could be difficult to achieve in present experiments. 

Before returning to the derivation of the master equation, we should re- 
mark that in realistic experimental situations, the squeezed modes cover only 
a small portion of the modes surrounding the atoms. The squeezing modes lie 
inside a cone of angle 6k < vr, and the modes outside the cone are in their ordi- 
nary vacuum state. In fact, the modes are in a finite temperature black-body 
state, which means that inside the cone the modes are in mixed squeezed vac- 
uum and black-body states. However, this is not a serious practical problem 
as experiments are usually performed at low temperatures where the black- 
body radiation is negligible. In principle, we can include the black-body 
radiation effect (thermal noise) to the problem replacing \D {uJk)f N (uk) in 
Eq. (|17|) by \D {ujk)f N (uk) + N, where N is proportional to the photon 
number in the black-body radiation. 

We now return to the derivation of the master equation for the density 
operator of the atomic system coupled to a squeezed vacuum field. First, we 
change the sum over ks into an integral 



E 



V 



^ poo p 

V / dukul / dQk ■ 



(20) 



Next, with the correlation functions ()17|) and after the rotating- wave ap- 
proximation (RWA) jHO], in which we ignore all terms oscillating at higher 
frequencies, 2uJi,uJi + Uj, the general master equation (fTB|) can be written as 

(t) = EZ=i { [SiX,, (t, r) , 5+] + [Sr, Xj, (t, r) S, 



+ 
+ 
+ 



Sj', yji (t, t) 



S+%it,T),Sr] + 
5+i:,,(t,r),5+] + [5+,i^,,(t,r)5; 
SrKl (t, r) , Sr] + [Sr , (t, r) Sj] } , (21) 



where the two-time operators are 

V 



Xi, (t, r) 



f dukule-'^^^-^^' f dnkY.Xij'^ {t,^ 

•' •' s=l 

(2vrc)3 / / E ^) 



(27rc)3 
V 



(22) 



s=l 
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with 



X 



V 



(27rc)= 



'i(2uJs—uJi~ujj)t 



\D{uk)VN{uk) + l 



X 



drp (t — r) e 



X 



xl"(«) 



/ drp {t — t) e 
Jo 



drp (t — t) e 



X ( 2 s — a; — o; J ) T 



(23) 



and is the sohd angle over which the squeezed vacuum field is propagated. 

The master equation (jTIj) with parameters and (^Hj) is quite general 
in terms of the matching of the squeezed modes to the vacuum modes and 
the bandwidth of the squeezed field relative to the atomic linewidths. The 
master equation is in the form of an integro-differential equation, and can 
be simplified by employing the Markov approximation In this approxi- 
mation the integral over the time delay r contains functions which decay to 
zero over a short correlation time Tc. This correlation time is of the order of 
the inverse bandwidth of the squeezed field, and the short correlation time 
approximation is formally equivalent to assume that squeezing bandwidths 
are much larger than the atomic linewidths. Over this short time-scale the 
density operator would hardly have changed from p (t) , thus we can replace 
p(t — t) hj p (t) in Eq. and extend the integral to infinity. Under these 
conditions, we can perform the integration over r and obtain [^0] 



lim / rfrp (t - r) e"" ^ p (t) 



tt6 (x) + i — 

X 



(24) 



where V indicates the principal value of the integral. Moreover, for squeez- 
ing bandwidth much larger than the atomic linewidths, we can approxi- 
mate the squeezing parameters and the mode function evaluated at Uk by 
their maximal values evaluated at Us, i.e., we can take N (uk) = N (us), 
M (tok) = M (us), and D {cuk) = D (w,). 

Finally, to carry out the polarization sums and integrals over dVLk in 
Eq. ()22|1 . we assume that the dipole moments of the atoms are parallel and 
use the spherical representation for the propagation vector k. The integral 
over dVLk contains integrals over the spherical angular coordinates 9 and (j). 
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The angle is formed by r^j and k directions, so we can write 



[sin 6 cos (f) , sin 6^ sin (j) , cos 6] 



(25) 



In this representation, the unit polarization vectors eg, and er^ may be chosen 



as 



[— COS 6 cos , — COS 6^ sin , sin 6] 
[sin , — cos 0,0], 



(26) 



and the orientation of the atomic dipole moments can be taken in the x 
direction 



Pi 



|/i,|[l ,0,0] 
1/2,1 [1,0,0] 



(27) 



With this choice of the polarization vectors and the orientation of the dipole 
moments, we obtain 



X,,(t,r) 
Yij {t, r) 

where 



(M) 



p{t) e 



-i{2L0s—L0i—UJj)t 



{2i 



with 



M(uJ^ 



N{uJs)\D{uj,)\'v{ 



1 - ^ (3 + cos^ 9s] cos9, 



(29) 



(30) 



and 9s is the angle over which the squeezed vacuum is propagated. 

The parameters Tij, which appear in Eq. (|^. are spontaneous emission 
rates, such that 



r = r 



(31) 



is the spontaneous emission rate of the ith atom, equal to the Einstein A 
coefficient for spontaneous emission, and 



r = r 

J- J- 



(^ ^ j) 



(32) 
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where 



3 

+ 



1 - (/^ ■ rijY 



sin (fconj) 
cos {k^rij) sin [korij] 



ihrij) {korijY 



(33) 



are collective spontaneous emission rates arising from the coupling between 
the atoms through the vacuum field [TTlll7lllHlinilE2jj. In the expression (j^ . 
fi = fii = Jij and fij are unit vectors along the atomic transition dipole 
moments and the vector -fij = fj — fi, respectively. Moreover, ko = uo/c, 
where Uq = + Ljj)/2, and we have assumed that — ujj) <^ luq. 

will 



The remaining parameters Q^f^ and ^^^^^ that appear in Eq. 



(M) 



contribute to the shifts of the atomic levels, and are given by 



(±) 



p 



rr 

J- 2 J- J 



27rix'o Jo u!k±ui 



-dujk 



(34) 



3 



and 



(M) 



P 



27ru;o -'o 



(35) 



where F (cj^ro / c) is given in Eq. with A;o replaced by uik / c, and 
replaced by ro = + r^. 

With the parameters ()28|). the master equation of the system of non- 
identical atoms in a broadband squeezed vacuum, written in the Schrodinger 
picture, reads 



dp 
dt 



1 



N 



- y 

2 ^ 



l + N(uJs 



pS^Sj + S^Sj- p — 2Sj pS^ 



1 JV 



AT 



*J=1 

N 



+^ E (r.. + ^^^) M (ujs) (pStSi + StSfp - 2S+pS, 



1 



where 



2N {ujs) + l] {^u^ - ^ 
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(-) 



(36) 



(37) 



represent a part of the intensity dependent Lamb shift of the atomic levels, 
while 



(-) 



^ j) 



(38) 



represents the vacuum induced coherent (dipole-dipole) interaction between 
the atoms. It is well known that to obtain a complete calculation of the Lamb 
shift, it is necessary to extend the calculations to a second-order multilevel 
Hamiltonian including electron mass renormalisation -63] ■ 

The parameters 6i are usually absorbed into the atomic frequencies Ui, 
by redefining the frequencies uji = uoi + 5i and are not often explicitly in- 
cluded in the master equations. The other parameters, fi-^^^ and do 
not appear as a shift of the atomic levels. One can show by the calculation 
of the integral appearing in Eq. (jH^j) that the parameter fij-^^ is negligibly 
small when the carrier frequency of the squeezed field is tuned close to the 
atomic frequencies IMj IHS 1^0] • On the other hand, the parameter VLij 
is independent of the squeezing parameters N {ojg) and M{ujs), and arises 
from the interaction between the atoms through the vacuum field. It can be 
seen that VLij plays a role of a coherent (dipole-dipole) coupling between the 
atoms. Thus, the collective interactions between the atoms give rise not only 
to the modified dissipative spontaneous emission but also lead to a coherent 
coupling between the atoms. 

Using the contours integration method, we find from Eq. (j38|) the explicit 
form of Vtij as [HI ill M., 67, 68] 



3 
4 

+ 



TiVj { - 1 - (/i ■ Tijf 



cos 



3 (yU ■ Ti 



sin {koTij) cos {k^r 



+ 



(39) 



and fijj, which both depend on the in- 



The collective parameters T 
teratomic separation, determine the collective properties of the multiatom 
system. In Fig. [H we plot Tij/^jTiTj and Vtij/ ^jTiTj as a function of r-jj/A, 
where A is the resonant wavelength. For large separations (r^ ^ A) the pa- 
rameters are very small {Tij = Qij ~ 0), and become important for r^j < A/2. 
For atomic separations much smaller than the resonant wavelength (the small 
sample model), the parameters attain their maximal values 



and 



3Jr,r, 



3 (fl ■ Tijf 



(40) 



(41) 



4 {kQTijf 

In this small sample model VLij corresponds to the quasistatic dipole-dipole 
interaction potential. 
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Figure 1: (a) Collective damping Tij/ JViVj and (b) the dipole-dipole inter- 



action Vtij / ^^TiT j as a function of rij/X for jj, ± Vij (solid line) and /i 
(dashed line). 



r 



1-3 



Equation ()36|) is the final form of the master equation that gives us an 
elegant description of the physics involved in the dynamics of interacting 
atoms. The collective parameters Tij and fijj, which arise from the mutual 
interaction between the atoms, significantly modify the master equation of 
a two-atom system. The parameter Tij introduces a coupling between the 
atoms through the vacuum field that the spontaneous emission from one of 
the atoms influences the spontaneous emission from the other. The dipole- 
dipole interaction term VLij introduces a coherent coupling between the atoms. 
Owing to the dipole-dipole interaction, the population is coherently trans- 
ferred back and forth from one atom to the other. Here, the dipole-dipole 
interaction parameter Vtij plays a role similar to that of the Rabi frequency 
in the atom-field interaction. 

For the next few sections, we restrict ourselves to the interaction of the 
atoms with the ordinary vacuum, M{ijJs) = N{uJs) = 0, and driven by an 
external coherent laser field. In this case, the master equation (j36j) can be 
written as 



dp i - 



1 



N 



2 E [pS^Sr + S^Srp - 2Sj pSt) , (42) 
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where 



N N 

Hs = nJ2{u, + 5^)S! + hJ2n,,S+Sr + HL, (43) 

i=l i^j 

and 

1 ^ 

Hl = --^E[^(^"'0'S'+e^(-^*+'^^)+H.c.] , (44) 

i=l 

is the interaction Hamiltonian of the atoms with a classical coherent laser 
field of the Rabi frequency (rj), the angular frequency ojl and phase 0^. 

Note that the Rabi frequencies of the driving field are evaluated at the 
positions of the atoms and are defined as |60| 

(^(r-;) = a = A?^-^Le*'^"7^, (45) 

where is the amplitude and is the wave vector of the driving field, 
respectively. The Rabi frequencies depend on the positions of the atoms and 
can be different for the atoms located at different points. For example, if the 
dipole moments of the atoms are parallel, the Rabi frequencies Vti and Vtj of 
two arbitrary atoms separated by a distance r^j are related by 

= ^].lMe^^L-^^- ^ (46) 

I A'-i I 

where Tij is the vector in the direction of the interatomic axis and \ fij\ = rij is 
the distance between the atoms. Thus, for two identical atoms (|/ij| = |/ij|), 
the Rabi frequencies differ by the phase factor exp{ikL ■ fij) arising from dif- 
ferent position coordinates of the atoms. However, the phase factor depends 
on the orientation of the interatomic axis in respect to the direction of prop- 
agation of the driving field, and therefore exp^iki ■ rij) can be equal to one, 
even for large interatomic separations rij. This happens when the direction 
of propagation of the driving field is perpendicular to the interatomic axis, 
ki ■ rij = 0. For directions different from perpendicular, ki ■ fij ^ 0, and then 
the atoms are in nonequivalent positions in the driving field, with different 
Rabi frequencies (fij 7^ fij). For a very special geometrical configuration of 
the atoms that are confined to a volume with linear dimensions that are much 
smaller compared to the laser wavelength, the phase factor exp(ifc2. 'rij) ~ 1, 
and then the Rabi frequencies are independent of the atomic positions. This 
specific configuration of the atoms is known as the small sample model or the 
Dicke model, and do not correspond in general to the experimentally realised 
atomic systems such as atomic beams or trapped atoms. 

The formalism presented here for the derivation of the master equation 
can be easily extended to the case of multi-level atoms [HHl EOl HH IZ2] and 
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atoms interacting with colour (frequency dependent) reservoirs |7^I7HI751 I7U] 
or photonic band-gap materials [771 EH]- Freedhoff has extended the 
master equation formahsm to electric quadrupole transitions in atoms. In 
the following sections, we will apply the master equations and to a 
wide variety of cases ranging from two identical as well as nonidentical atoms 
interacting with the ordinary vacuum to atoms driven by a laser field and 
finally to atoms interacting with a squeezed vacuum field. 

2.2 Quantum jump approach 

The master equation is a very powerful tool for calculations of the dynamics 
of Markovian systems which assume that the bandwidth of the vacuum field 
is broadband. The Markovian master equation leads to linear differential 
equations for the density matrix elements that can be solved numerically or 
analytically by the direct integration. 

An alternative to the master equation technique is quantum jump ap- 
proach. This technique is based on quantum trajectories [50^ that are equiv- 
alent to the Monte Carlo wave-function approach 13^ . and has been 
developed largely in connection with problems involving prediction of all 
possible evolution trajectories of a given system. This approach can be used 
to predict all evolution trajectories of a single quantum system which stochas- 
tically emits photons. Our review of this approach will concentrate on the 
example considered by Beige and Hegerfeldt [80 of two identical two-level 
atoms interacting with the three-dimensional EM field whose the modes are 
in the ordinary vacuum states. 

In the quantum jump approach it is assumed that the probability density 
for a photon emission is known for all times t, and therefore the state of 
the atoms changes abruptly. After one photon emission the system jumps 
into another state, which can be determined with the help of the so called 
reset operator. The continuous time evolution of the system between two 
successive photon emissions is determined by the conditional Hamiltonian H^. 
Suppose that at time to the state of the combined system of the atoms and 
EM field is given by 



where p is the density operator of the atoms and |0) is the vacuum state of 
the field. After a time At a photon is detected and then the state of the 
system changes to 



where P = 1 — |0) (0| is the projection onto the one photon space, and 



^)(^l = |0)p(0| , 



(47) 



VUi (to + At, to) |0) p (0| Uj (to + At, to) V 



(48) 



?7,(t,to)=e-^^W(*-*°) 



(49) 
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is the evolution operator with the Hamiltonian V{t) given in Eq. (jH}. 

The non-normahsed state of the atomic system, denoted as R{p)At, is 
obtained by taking trace of Eq. PHj) over the field states 



(50) 



R{p)At = Tip (vUi (to + At, to) |0) p (0| U\ (to + At, to) V 



where R{p) is called the non-normalised reset state and the corresponding 
operator R is called the reset operator. 

Using the perturbation theory and Eq. (jH)), we find the explicit form of 
R{p) for the two-atom system as 



\ (Ci*2 + C21) s^psi + \ {c,2 + c;,) s.pst 



+T[S^pSt + S,pS^ 



(51) 



where 



Ci' 



2 



{p-r. 



+ 



3 (/i ■ rijf 



koVij 



(52) 



Note that ReCjj = Tij and ImCij = 2Qij, where and Qij are the collective 
atomic parameters, given in Eqs. and (jSHl), respectively. 

The time evolution of the system under the condition that no photon is 
emitted is described by the conditional Hamiltonian He, which is found from 
the relation 



1 - -HrAt 



(0|[/,(to + At,to)|0) 



(53) 



where At is a short evolution time such that At < l/F. Using second order 
perturbation theory, we find from Eq. ()53|) that the conditional Hamiltonian 
for the two-atom system is of the form 



2i 



r [Si Si + ^2 5*2 ) + C12S1 S2 + C2l5'2 Si 



(54) 



Hence, between photon emissions the time evolution of the system is given 
by an operator 



Uc (to + At, to 



(55) 



which is nonunitary since He is non-Hermitian, and the state vector of the 
system is 



I* 



At/ 



f/e(to + At,to)|^o) • 



(56) 
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Then, the probabihty to detect no photon until time t is given by 

P{t;\^o)) = |?7c(t,to)|^o)f . (57) 

The probabihty density Wi (t; I^E'o)) of detecting a photon at time t is defined 
as 

«;i(t;|^o)) = -^P(t;|^o)), (58) 

and is often caUed the waiting time distribution. 

The resuhs fjHTjl and show that in the quantum jump method one 
calculates the times of the photon detection stochastically. Starting at t = to 
with a pure state, the state develops according to JJc until the first emission at 
some time ti, determined from the waiting time Wi. Then the state is reset, 
according to Eq. ()5H1 , to a new density matrix and the system evolves again 
according to Uc until the second emission appearing at some time ^2, and the 
procedure repeats until the final time tn- In this way, we obtain a set of tra- 
jectories of the atomic evolution. The ensemble of such trajectories yields to 
equations of motion which are solved using the standard analytical or numer- 
ical methods. As a practical matter, individual trajectories are generally not 
observed. The ensemble average over all possible trajectories leads to equa- 
tions of motion which are equivalent to the equations of motion derived from 
the master equation of the system. Thus, the quantum jump approach is 
consistent with the master equation method. However, the advantage of the 
quantum jump approach over the master equation method is that it allows to 
predict all possible trajectories of a single system. Using this approach, it has 
been demonstrated that environment induced measurements can assist in the 
realization of universal gates for quantum computing jTH|. Cabrillo et al. ^] 
have applied the method to demonstrate entangling between distant atoms 
by interference. Schon and Beige have demonstrated the advantage of 
the method in the analysis of a two- atom double-slit experiment. 



3 Entangled atomic states 

The modification of spontaneous emission by the collective damping and 
in particular the presence of the dipole-dipole interaction between the atoms 
suggest that the bare atomic states are no longer the eigenstates of the atomic 
system. We will illustrate this on a system of two identical as well as non- 
identical atoms, and present a general formalism for diagonalization of the 
Hamiltonian of the atoms in respect to the dipole-dipole interaction. 

In the absence of the dipole-dipole interaction and the driving laser field, 
the space of the two-atom system is spanned by four product states 

1^1)1^2) , |ei)|^2) , 1^1) |e2) , |ei)|e2) , (59) 
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with corresponding energies 



^gg = ~^^o , = -hA , Eg^ = hA , E^e = huo , (60) 

where uq = ^ {ui + UJ2) and A = ^ {uj2 — uJi). 

The product states \ei)\g2) and \gi)\e2) form a pair of nearly degenerated 
states. When we include the dipole-dipole interaction between the atoms, 
the product states combine into two linear superpositions (entangled states), 
with their energies shifted from ±^A by the dipole-dipole interaction energy. 
To see this, we begin with the Hamiltonian of two atoms including the dipole- 
dipole interaction 

Haa = E huJiS^ + hJ2 ^^JStS- . (61) 
i=l i^j 

In the basis of the product states ()59|) . the Hamiltonian ()6ip can be written 
in a matrix form as 



Han = h 



f -uo \ 

-A ^]l2 

fii2 A 

V ujo / 



(62) 



Evidently, in the presence of the dipole-dipole interaction the matrix ()62p 
is not diagonal, which indicates that the product states are not the 
eigenstates of the two-atom system. We will diagonalize the matrix (j62j) 
separately for the case of identical (A = 0) and nonidentical (A 7^ 0) atoms 
to find eigenstates of the systems and their energies. 

3.1 Entangled states of two identical atoms 

Consider first a system of two identical atoms (A = 0). In order to find 
energies and corresponding eigenstates of the system, we have to diagonalize 
the matrix ()62p . The resulting energies and corresponding eigenstates of the 
system are ^Ul HZj 

1^) = 1^1)1^2) , 

\s) = ^{\ei)\g2) + \gi)\e2)) , 

|a) = ^ (|ei)|c/2) - I^i)|e2)) , 
|e) = |ei)|e2) . (63) 

The eigenstates ()63p. first introduced by Dicke JU]; are known as the 
collective states of two interacting atoms. The ground state \g) and the upper 





= -'hujQ , 






Ea 


= -TiVti2 


Ee 
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Figure 2: Collective states of two identical atoms. The energies of the sym- 
metric and antisymmetric states are shifted by the dipole-dipole interaction 
Qi2. The arrows indicate possible one-photon transitions. 

state |e) are not affected by the dipole-dipole interaction, whereas the states 
|s) and I a) are shifted from their unperturbed energies by the amount ±f2i2, 
the dipole-dipole energy. The most important property of the collective states 
|s) and I a) is that they are an example of maximally entangled states of the 
two-atom system. The states are linear superpositions of the product states 
which cannot be separated into product states of the individual atoms. 

We show the collective states of two identical atoms in Fig. |21 It is seen 
that in the collective states representation, the two-atom system behaves as 
a single four- level system, with the ground state \g), the upper state |e), 
and two intermediate states: the symmetric state \s) and the antisymmetric 
state I a). The energies of the intermediate states depend on the dipole- 
dipole interaction and these states suffer a large shift when the interatomic 
separation is small. There are two transition channels |e) \s) \g) and 
|e) — i> \a) \g), each with two cascade nondegenerate transitions. For two 
identical atoms, these two channels are uncorrelated, but the transitions in 
these channels are damped with significantly different rates. To illustrate 
these features, we transform the master equation ()42|) into the basis of the 
collective states (jUHj) . We define collective operators Aij = where 
i,j = e, a, s, g, that represent the energies (z = j) of the collective states and 
coherences (z 7^ j). Using Eq. (jHS)), we find that the collective operators are 
related to the atomic operators Sf through the following identities 




1 



(A 



V2 



■es 



■ea ~l~ -^sg ~l~ -^ag) i 
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71 ^^^^ + 



(64) 



Substituting the transformation identities into Eq. (j42j), we find that in 
the basis of the collective states the master equation of the system can be 
written as 



Ft' 



I r 



n 



Has, P 



(65) 



where 



Has = h [C^O {Aee " Agg) + Vti2 {A^s - Aaa)] 

- ^ (^1 + ^2) [{Aes + Asg) e^(-^*+<^^) + H.c. 

^ (^2 - ^^l) \{Aea - Aag) e^(-^*+^-) + H.C. 



is the Hamiltonian of the interacting atoms and the driving laser field, 



i (r + Tu) {{Aee + Ass) P + P {Aee + As 



(66) 




~ 2 {Age + Ags) p {Aes + Asg)} 



(67) 



describes dissipation through the cascade |e) 
the symmetric state \s), and 



|s) —>■ \g) channel involving 




1 



(r - Tia) {{Aee + Aaa) P + P {Aee + Aa 



~ 2 {Aae — Aga) f) {Aea ~ ^ag)} 



(68) 



describes dissipation through the cascade |e) \a) \g) channel involving 
the antisymmetric state \a). 

We will call the two cascade channels |e) — * \s) \g) and |e) — > \a) — > l^f) 
as symmetric and antisymmetric transitions, respectively. The first term in 
Has is the energy of the collective states, while the second and third terms 
are the interactions of the laser field with the symmetric and antisymmetric 
transitions, respectively. One can see from Eqs. ()65|) - (jUHj) that the symmet- 
ric and antisymmetric transitions are uncorrelated and decay with different 
rates; the symmetric transitions decay with an enhanced (superradiant) rate 
(r + ri2), whereas the antisymmetric transitions decay with a reduced (sub- 
radiant) rate (F — ri2). For F = F12, which appears when the interatomic 
separation is much smaller than the resonant wavelength, the antisymmet- 
ric transitions decouple from the driving field and does not decay. In this 
case, the antisymmetric state is completely decoupled from the remaining 
states and the system decays only through the symmetric channel. Hence, 
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for = r the system reduces to a three-level cascade system, referred to 
as the small-sample model or two-atom Dicke model ^21 E3 HH] • The model 
assumes that the atoms are close enough that we can ignore any effects result- 
ing from different spatial positions of the atoms. In other words, the phase 
factors exp{ik ■•f'i) are assumed to have the same value for all the atoms, and 
are set equal to one. This assumption may prove difficult in experimental 
realization as the present atom trapping and cooling techniques can trap two 
atoms at distances of the order of a resonant wavelength IT^ IT^ HE] . 
At these distances the collective damping parameter differs significantly 
from r (see Fig. Q), and we cannot ignore the transitions to and from the 
antisymmetric state. We can, however, employ the Dicke model to spatially 
extended atomic systems. This could be achieved assuming that the obser- 
vation time of the atomic dynamics is shorter than F"^. The antisymmetric 
state \a) decays on a time scale ~ (F — Fi2)~^, which for F12 ~ F is much 
longer than F~^. On the other hand, the symmetric state decays on a time 
scale ~ (F -|- Fi2)~^, which is shorter than F~^. Clearly, if we consider short 
observation times, the antisymmetric state does not participate in the dy- 
namics and the system can be considered as evolving only between the Dicke 
states. 

Although the symmetric and antisymmetric transitions of the collective 
system are uncorrelated, the dynamics of the four-level system may be signif- 
icantly different from the three-level Dicke model. As an example, consider 
the total intensity of the fluorescence field emitted from a two-atom system 
driven by a resonant coherent laser field {ujl = ujq). We make two simplifying 
assumptions in order to obtain a simple analytical solution: Firstly, we limit 
our calculations to the steady-state intensity. Secondly, we take ■ ri2 = 
that corresponds to the direction of propagation of the driving field perpen- 
dicular to the interatomic axis. We emphasize that these assumptions do not 
limit qualitatively the physics of the system, as experiments are usually per- 
formed in the steady-state, and with ki ■ fi2 = the interatomic separation 
ri2 may still be any size relative to the resonant wavelength. 

We consider the radiation intensity I{R,t) detected at a point R at the 
moment of time t. If the detection point R is in the far-field zone of the 
radiation emitted by the atomic system, then the intensity can be expressed 
in terms of the first-order correlation functions of the atomic dipole operators 

as Ha Hi 

2 

/ (i?, t) = u{R) {St (t - R/c) SJ (t - R/c))e'^^-''^^ , (69) 

where 

u{R) = (tu^/iV2i?^cVeo) sin^ cp (70) 

is a constant which depends on the geometry of the system, ip the angle 
between the observation direction R = RR and the atomic dipole moment 
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/i. 

On integrating over all directions, Eq. ()69|) yields the total radiation in- 
tensity given in photons per second as 

/ W = E ^^ASt it - R/c) S- it - R/c)) . (71) 

The atomic correlation functions, appearing in Eq. (IHJ, are found from the 
master equation (jl^ . There are, however, two different steady-state solu- 
tions of the master equation ()42|) depending on whether the collective damp- 
ing rates Tis = T or Tis ^ T [83, 84. 

For 7^ r and ■ ri2 = 0, the steady-state solutions for the atomic 
correlation functions are 

AD ' 

StS,) = {S^S^) = (72) 



where 



D = n^ + [n^ + nl,) + Ir^ (r + r^f . (73) 

If we take = F and = 0, that corresponds to the two-atom Dicke 
model, the steady-state solutions for the atomic correlation functions are of 
the following form 



2D' 



StS,) = {S^S^) = — , (74) 

where 

D' = 3Q'^ + 4F2fi2 + . (75) 

In the limit of a strong driving field, F, the steady-state total radia- 
tion intensity from the two-atom Dicke model is equal to 4F/3. However, for 
the spatially separated atoms Igs = limt^oo I [t) = F, which is twice of the 
intensity from a single atom jHS] ■ There is no additional enhancement of the 
intensity. 

Note that in the limit of ri2 0, the steady-state solution (f7^ does not 
reduce to that of the Dicke model, given in Eq. (ff^ . This fact is connected 
with conservation of the total spin S"^, that S"^ is a constant of motion for the 
Dicke model and 5*^ not being a constant of motion for a spatially extended 
system of atoms jHSl IHl] ■ We can explain it by expressing the square of the 
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total spin of the two-atom system in terms of the density matrix elements of 
the collective system as 

(t) = 2 - 2paa (t) . (76) 

It is clear from Eq. ()76p that S"^ is conserved only in the Dicke model, in 
which the antisymmetric state is ignored. For a spatially extended system 
the antisymmetric state participates fully in the dynamics and S"^ is not 
conserved. The Dicke model reaches steady state between the triplet states 
|e), and \g), while the spatially extended two-atom system reaches steady 
state between the triplet and the antisymmetric states. 

Amin and Cordes |H7] calculated the total radiation intensity from an A^- 
atom Dicke model and showed the intensity is A^(A^ -|- 2)/3 times that for a 
single atom, which they called " scaling factor" . The above calculations show 
that the scaling factor is characteristic of the small sample model and does 
not exist in spatially extended atomic systems. Thus, in physical systems the 
antisymmetric state plays important role and as we have shown its presence 
affects the steady-state fluorescence intensity. The antisymmetric state can 
also affect other phenomena, for example, photon antibunching jHH], and 
purity of two-photon entangled states, that is discussed in Sec. M 

3.2 Collective states of two nonidentical atoms 

For two identical atoms, the dipole-dipole interaction leads to the maxi- 
mally entangled symmetric and antisymmetric states that decay indepen- 
dently with different damping rates. Furthermore, in the case of the small 
sample model of two atoms the antisymmetric state decouples from the ex- 
ternal coherent field and the environment, and consequently does not decay. 
The decoupling of the antisymmetric state from the coherent field prevents 
the state from the external coherent interactions. This is not, however, an 
useful property from the point of view of quantum computation where it is 
required to prepare entangled states which are decoupled from the external 
environment and simultaneously should be accessible by coherent processes. 
This requirement can be achieved if the atoms are not identical, and we will 
discuss here some consequences of the fact that the atoms could have differ- 
ent transition frequencies or different spontaneous emission rates. To make 
our discussion more transparent, we will concentrate on two specific cases: 
(1) A ^ and Fi = Fa, and (2) A = and Fi ^ Fa. 

3.2.1 The case A ^ and Fi = Fs 

When the atoms are nonidentical with different transition frequencies, the 
states (|63p are no longer the eigenstates of the Hamiltonian (|6U|) . The diago- 
nalization of the matrix (jU^ with A 7^ leads to the following energies and 
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corresponding eigenstates jHS] 



Eg = -hUJO , l^f) = |fi'l)|fi'2) , 

= hw , \s') = P\ei)\g2) + a|fl'i)|e2) , 

Ea' = -hw , \a') = a\ei)\g2) - P\gi)\e2) , 

Ee = huo , |e) = |ei)|e2) , (77) 



where 



a 



, p= '\ , w = ^ni2 + ^', (78) 



and ci = A + ^^l^ + A2. 

The energy level structure of the collective system of two nonidentical 
atoms is similar to that of the identical atoms, with the ground state \g), 
the upper state |e), and two intermediate states and \a'). The effect 
of the frequency difference A on the collective atomic states is to increase 
the splitting between the intermediate levels, which now is equal io w = 



Vl\2 + A^. However, the most dramatic effect of the detuning A is on the 
degree of entanglement of the intermediate states |s') and \a') that in the case 
of nonidentical atoms the states are no longer maximally entangled states. 
For A = the states |s') and \a') reduce to the maximally entangled states 
\s) and I a), whereas for A ^ VL12 the entangled states reduce to the product 
states \ei)\g2) and —\gi)\e2), respectively. 

Using the same procedure as for the case of identical atoms, we rewrite 
the master equation in terms of the collective operators Aij = \i){j\, 
where now the collective states \i) are given in Eq. (j77|) . First, we find that 
in the case of nonidentical atoms the atomic dipole operators can be written 
in terms of the linear combinations of the collective operators as 

= (3Aes' + aAea' + aAs>g - (3Aa'g . (79) 

Hence, in terms of the collective operators Aij, the master equation takes the 
form 

d i r - 1 

Hs',p +Cp, (80) 



where 



Hs' = h[LJo{Aee- Agg) +w{As's' - Aa'a')] 

- I { [(afil + ^^2) Aes' + + 01^2) A,g\ e*(-^*+<^^) 

+ [(«fi2-/3fil)Aa'-(/5f^2-«f^l)^age^^-^*+*^)+H.C.} (81) 
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is the Hamiltonian of the system in the collective states basis, and the Liou- 
ville operator Cp describes the dissipative part of the evolution. The dissi- 
pative part is composed of three terms 



where 



= -r^' { {Aee + As's' )p + p{Aee + A^'s' ) 

-2 {As'epAes' + Ags'pAs'g)} 

- {a(3T + Tu) {As'epAs'g + Ags'pAes') , (83) 

■^Pj = -ra' { {Aee + Aa'a' )P + P (^ee + Aa'a' ) 
-2 {Aa'epAea' + Aga'pAa'g)] 

-{a(3T -V 12) {Aa'epAa'g + Aga'pAea') , (84) 



and 



i^Pj = -Ta's' {{Aa's' + As'a') P + pica's' + As'a') 

- 2 {Aga'pAs'g + Ags'pAa'g + As'epAea' + Aa'epAes')} 
+ - fi^) V {Aa'epAs'g + Ags'pA^a' 

+As'epAa'g + Aga'pAes'} , (85) 

with the damping coefficients 

Ts' = ^ (r + 2a/3ri2) , r,, = ^ (r - 2a(3ru) , 

Tav = ^{a'-f3')ru . (86) 

The dissipative part of the master equation is very extensive and unlike 
the case of identical atoms, contains the interference term between the sym- 
metric and antisymmetric transitions. The terms (jHHjl and ()84|) describes 
spontaneous transitions in the symmetric and antisymmetric channels, re- 
spectively. The coefficients F,./, and r^/ are the spontaneous emission rates 
of the transitions. The interference term (j85p results from spontaneously 
induced coherences between the symmetric and antisymmetric transitions. 
This term appears only in systems of atoms with different transition frequen- 
cies (A 7^ 0), and reflects the fact that, as the system decays from the state 
|s'), it drives the antisymmetric state, and vice versa. Thus, in contrast to 
the case of identical atoms, the symmetric and antisymmetric transitions are 
no longer independent and are correlated due to the presence of the detuning 
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Figure 3: The spontaneous emission damping rate Ta' as a function of A 
for /2 ± fi2, and different interatomic separations: r^/X = 0.05 (solid line), 
ru/X = 0.1 (dashed line), ri2/A = 0.5 (dashed-dotted line). 

A. Moreover, for nonidentical atoms the damping rate of the antisymmetric 
state cannot be reduced to zero. In the case of interatomic separations much 
smaller than the optical wavelength (the small sample model), the damping 
rate reduces to 

Ta' = ^T{a-pf , (87) 

that is different from zero, unless A = 0. 

In Fig. 121 we plot the damping rate Ta' as a function of A for different 
interatomic separations. The damping rate vanishes for A = independent 
of the interatomic separation, but for small interatomic separations there is 
a significant range of A for which F^/ <C F. 

3.2.2 The case A = and Fi ^ Fa 

The choice of the collective states (f77|) as a basis leads to a complicated 
dissipative part of the master equation. A different choice of collective states 
is proposed here, which allows to obtain a simple master equation of the 
system with only the uncorrelated dissipative parts of the symmetric and 
antisymmetric transitions ^H]- Moreover, we will show that it is possible to 
create an entangled state in the system of two nonidentical atoms which can 
be decoupled from the external environment and, at the same time, the state 
exhibits a strong coherent coupling with the remaining states. 
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To illustrate this, we introduce superposition operator Sf" and which 
are linear combinations of the atomic dipole operators Sf" and as 

Sj' = uSi + , = u*Si + v*S2 , 

= vSt-uS+, S- = v*S^-u*S2 , (88) 

where u and v are the transformation coefficients which are in general com- 
plex numbers. The coefficients satisfy the condition 

\uf + \vf = l. (89) 

The operators 5^ and 5*^ represent, respectively, symmetric and anti- 
symmetric superpositions of the atomic dipole operators. In terms of the 
superposition operators, the dissipative part of the master equation (jl^ can 
be written as 

Cp = -Tss {S+S;p + pS+S- - 2S;pS+) 
— Taa (^S'^S^ p + pS^S^ — 2S^ pS'^ 
— Tsa (^S^S^ p + pS'^S^ — 2S^ pS^ 

-r,, (s:s;p + pS^S; - 2S;pS: ) , (90) 



where the coefficients F^n are 

Tss = \ufTi + \vfr2 + {uv* + u*v)ri2, 

Taa = \vfTi + \ufT2-{uV* + U*v)Ti2, 
Tas = UV*ri - U* VT 2 - (juf - {vf'jT 12, 

Tsa = U*vTi - UV*T 2 - {\uf ~ Ivf'jT 12. (91) 

The first two terms in Eq. ()9()j) are familiar spontaneous emission terms of 
the symmetric and antisymmetric transitions, and the parameters Tgg and 
Taa are spontaneous emission rates of the transitions, respectively. The last 
two terms are due to coherence between the superposition states and the 
parameters Tas and Tsa describes cross-damping rates between the superpo- 
sitions. 

If we make the identification 



then the damping coefficients ()9H) simplify to 

1 , , (ri2 - VT\r: 

r,, = -(ri + r2)+ ^ 



r 

aa 



2' ' ^' Ti + Ta 

i(ri-r2) (v/rj^-ri2 



2 Ti 
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(93) 



When the damping rates of the atoms are equal (Fi = r2), the cross-damping 
terms Tas and vanish. Furthermore, if = ^/TiT^ then the spontaneous 
emission rates Taa, ^as and Tsa vanish regardless of the ratio between the Fi 
and r2. In this case, which corresponds to interatomic separations much 
smaller than the optical wavelength, the antisymmetric superposition does 
not decay and also decouples from the symmetric superposition. 

An interesting question arises as to whether the nondecaying antisymmet- 
ric superposition can still be coupled to the symmetric superposition through 
the coherent interactions Q12 and fl contained in the Hamiltonian Hg- These 
interactions can coherently transfer population between the superpositions. 
To check it, we first transform the Hamiltonian PHj) into the interaction pic- 
ture and next rewrite the transformed Hamiltonian in terms of the Sf and 

operators as 

= -hAL[{StS; + S:S-) + iv*u-vu*){S^S--S:S;)' 

+nn,2 {{vu* + v*u) {sts; - s^s-) 
+ [\v\'~\u\') {sts- + s:s;)} 

~h[{uni + vn2)st + {vni-un2)s^ + 11.C.] , (94) 

where Al = ujl ~ t^o- 

In the above equation, the first term arises from the atomic Hamiltonian 
and shows that in the absence of the interatomic interactions the symmetric 
and antisymmetric states have the same energy. The second term in Eq. ()94|). 
proportional to the dipole-dipole interaction between the atoms, has two 
effects on the dynamics of the symmetric and antisymmetric superpositions. 
The first is a shift of the energies and the second is the coherent interaction 
between the superpositions. It is seen from Eq. (jU^ that the contribution 
of Q12 to the coherent interaction between the superpositions vanishes for 
Fi = F2 and then the effect of Qu is only the shift of the energies from their 
unperturbed values. Note that the dipole-dipole interaction Q12 shifts the 
energies in the opposite directions. The third term in Eq. (j94j) represents 
the interaction of the superpositions with the driving laser field. We see 
that the symmetric superposition couples to the laser field with an effective 
Rabi frequency proportional to uQi + VQ2, whereas the Rabi frequency of 
the antisymmetric superposition is proportional to vQi — UQ2 and vanishes 
for vQi = ufl2- 

Alternatively, we may write the Hamiltonian (|94j) in a more transparent 
form which shows explicitly the presence of the coherent coupling between 
the symmetric and antisymmetric states 



H. = -fi 

-ifi 
2 



{Al - A') StS; + {Al + A') S:S- + A^S^S- + AIS^S, 

{uni + vn2)st + {vni-un2)s+ + R.C] , (95) 
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where A' and Ac are given by 

A' = {vu* + v*u)Qi2 , = (\u\^ - Ivl"^) Qi2 + {v*u - vu*) Al . (96) 

The parameters A' and A^ allow us to gain physical insight into how the 
dipole-dipole interaction Q12 and the unequal damping rates Fi 7^ r2 can 
modify the dynamics of the two-atom system. The parameter A' appears 
as a shift of the energies of the superposition systems, while Ac determines 
the magnitude of the coherent interaction between the superpositions. For 
identical atoms the shift A' reduces to Q12 that is the dipole-dipole interaction 
shift of the energy levels. In contrast to the shift A', which is different from 
zero for identical as well as nonidentical atoms, the coherent coupling Ac can 
be different from zero only for nonidentical atoms. 

Thus, the condition = a/TiT^ for suppression of spontaneous emission 
from the antisymmetric state is valid for identical as well as non-identical 
atoms, whereas the coherent interaction between the superpositions appears 
only for nonidentical atoms with different spontaneous damping rates. 

It should be noted that this treatment is valid with only a minor mod- 
ification for a number of other schemes of two-atom systems. For example, 
it can be applied to the case of two identical atoms that experience different 
intensities and phases of the driving field [Hill EB 102] • 

In what follows, we will illustrate how the interference term in the master 
equation of two nonidentical atoms results in quantum beats and transfers 
of the population to the antisymmetric state even if the antisymmetric state 
does not decay. Of particular interest is the temporal dependence of the total 
radiation intensity of the fluorescence field emitted by two interacting atoms. 

4 Quantum beats 

The objective of this section is to give an account of interference effects result- 
ing from the direct correlations between the symmetric and antisymmetric 
states. We will first analyse the simplest model of spontaneous emission 
from two nonidentical atoms and consider the time dependence of the total 
radiation intensity. After this, we will consider the time evolution of the 
fluorescence intensity emitted by two identical atoms that are not in the 
equivalent positions in the driving field. 

4.1 Quantum beats in spontaneous emission from two 
nonidentical atoms 

For two nonidentical atoms the master equation ()42|) . in the absence of the 
driving field {Qi = 0), leads to a closed set of five equations of motion for the 
expectation values of the atomic dipole operators [8R . This set of equations 
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can be written in a matrix form as 



AX {t) 



where X (t) is a column vector with components 

Xi = {St{t)S^{t)) , X, = {SUt)S^{t)) , 
X, = {Snt)S,it)), X,= {S^it)S^{t)) , 

X, = {snt)snt)s^{t)s,{t)) , 

and A is the 5x5 matrix 
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(97) 



(98) 



(99) 



with n = —\{Ti2 + ^^12) and = Fi + r2. 

It is seen from Eq. (|99p that the equation of motion for the second-order 
correlation function {Si (t) S2 (t) Si (t) S2 (t)) is decoupled from the remaining 
four equations. This allows for an exact solution of the set of equations ()96|) . 
The exact solution is given in Ref. jHn]. Here, we will focus on two special 
cases of A 7^ 0, Fi = r2 and A = 0, Fi 7^ r2, and calculate the time evolution 
of the total fluorescence intensity, defined in Eq. (f7T|) . We will assume that 



initially {t = 0) atom 
its ground state \g2)- 



'1" was in its excited state |ei) and atom "2" was in 



4.1.1 The case A ^ 0, Fi = r2 = F and ^12 > A 

In this case the atoms have the same spontaneous damping rates but different 
transition frequencies that, for simplicity, are taken much smaller than the 
dipole-dipole interaction potential. In this limit, the approximate solution of 
Eq. (P7j) leads to the following total radiation intensity 



I it) 



-rt 



A 



2n 



-F12 cos 2wt + F cosh ri2t — F12 sinh ri2t 



12 



(100) 



where w = y^i2 + 

The total radiation intensity exhibits sinusoidal modulation (beats) su- 
perimposed on exponential decay with the damping rates F ± ri2. The 
amplitude of the oscillations is proportional to A and vanishes for identical 
atoms. The damping rate F + ri2 describes the spontaneous decay from 
the state \s') to the ground state \g), while F — ri2 is the decay rate of the 
I a') — > \g) transition. The frequency 2w of the oscillations is equal to the 
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frequency difference between the \s') and \a') states. The oscillations reflect 
the spontaneously induced correlations between the |s') \g) and \a') — > \g) 
transitions. According to Eq. ()86|) the amplitude of the spontaneously in- 
duced correlations is equal to Ta's', which in the limit of fli2 3> A reduces 
to Ta's' = Ari2/(2i7i2). Hence, the amplitude of the oscillations appearing 
in Eq. ()1UU|) is exactly equal to the amplitude of the spontaneously induced 
correlations. Fig. |3] shows the temporal dependence of the total radiation 
intensity for interatomic separation ri2 = A/12, Fi = F2, /i -L r, and different 
A. As predicted by Eq. (jlUUp . the intensity exhibits quantum beats whose 
the amplitude increases with increasing A. Moreover, at short times, the 
intensity can become greater than its initial value I (0). This effect is known 
as a superradiant behavior and is absent in the case of two identical atoms. 
Thus, the spontaneously induced correlations between the \s') \g) and 
\a') \g) transitions can induce quantum beats and superradiant effect in 
the intensity of the emitted field. 

The superradiant effect is characteristic of a large number of atoms jHSHHH 
1^ , and it is quite surprising to obtain this effect in the system of two atoms. 
Coffey and Friedberg (HH] and Richter have shown that the superradiant 
effect can be observed in some special cases of the atomic configuration of a 
three-atom system. Blank et al. [HH] have shown that this effect, for atoms lo- 
cated in an equidistant linear chain, appears for at least six atoms. Recently, 
DeAngelis et al. have experimentally observed the superradiant effect in 
the radiation from two identical dipoles located inside a planar symmetrical 
microcavity. 

Quantum beats predicted here for spontaneous emission from two non- 
identical atoms are fully equivalent to the quantum beats predicted recently 
by Zhou and Swain |lUUj in a single three-level V system with correlated 
spontaneous transitions. For the initial conditions used here that initially 
only one of the atoms was excited, the initial population distributes equally 
between the states \s') and \a'). Since the transitions are correlated through 
the dissipative term Ta's'-, the system of two nonidentical atoms behaves as 
a three-level V system with spontaneously correlated transitions. 

4.1.2 The case of A = 0, Fi ^ F2 and 9.12 > Fi, F2 

We now wish to show how quantum beats can be obtained in two nonidenti- 
cal atoms that have equal frequencies but different damping rates. According 
to Eqs. (|93|) and (|96p . the symmetric and antisymmetric transitions are cor- 
related not only through the spontaneously induced coherences F^g, but also 
through the coherent coupling Ac. One can see from Eq. (j^Hjl that for small 
interatomic separations Tas ~ 0. However, the coherent coupling parameter 
Ac, which is proportional to i7i2, is very large, and we will show that the 
coherent coupling Ac can also lead to quantum beats and the superradiant 
effect. In the case of A = 0,Fi 7^ F2 and ^12 ^ Fi,F2, the approximate 
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Figure 4: Time evolution of the total radiation intensity for = A/12, Fi = 
F2,/i -L f, and different A: A = (solid line), A = — 2F (dashed line), 
A = — 3F (dashed-dotted line). 

solution of Eq. ^7\i leads to the following expression for the total radiation 
intensity 

Ht) = e-^(r^+r^)*|^(Fi-F2)cos21]i2t 

+ ^(Fi + F2) cosh Fi2t-Fi2 sinh Fi2t| . (101) 

The intensity displays quantum-beat oscillations at frequency 2f2i2 corre- 
sponding to the frequency splitting between the |s') and \a') states. The 
amplitude of the oscillations is equal to (Fi — F2) /2 that is proportional to 
the coherent coupling A^. For Fi = F2 the coherent coupling parameter 
Ac = and no quantum beats occur. In this case the intensity exhibits pure 
exponential decay. This is shown in Fig. where we plot the time evolution 
of I (t) for interatomic separation = A/12, and different ratios F2/F1. 
Similar to the case discussed in Sec. 14.1.11 the intensity exhibits quantum 
beats and the superradiant effect. For = A/12 the collective damping 
F12 ~ v^FiF^, and then the parameter F^s ~ 0, indicating that the quantum 
beats and the superradiant effect result from the coherent coupling between 
the \s') and \a') states. 
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Figure 5: Time evolution of the total radiation intensity for r 12 = A/12, A = 
0,/i ± fi2, and different r2/ri: r2/ri = 1 (solid line), r2/ri = 2.5 (dashed 
line), r2/ri = 5 (dashed-dotted line). 

4.1.3 Two identical atoms in nonequivalent positions in a driving 
field 

Quantum beats and superradiant effect induced by interference between dif- 
ferent transitions in the system of two nonidentical atoms also occur in other 
situations. For example, quantum beats can appear in a system of two iden- 
tical atoms that experience different amplitude or phase of a coherent driving 

field inniEi]. 

Consider the Hamiltonian pij) of the interaction between coherent laser 
field and two identical atoms. In the interaction picture, the Hamiltonian 
can be written as 

Hl = -h[{niS+ + n2S^) + H.C.] , (102) 

where Qi is the Rabi frequency of the driving field at the position of the ith 
atom. 

For the atoms in a running- wave laser field with ki ■ fi ^ 0, the Rabi 
frequency is a complex parameter, which may be written as 

fii = ^(^^^-^^ , (103) 

where f2 = ■ Ei\/fi is the maximum Rabi frequency and is the wave 
vector of the driving field. Thus, in the running-wave laser field the atoms 
experience different phases of the driving field. 
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For the atoms in a standing-wave laser field and ki ■ ri ^ 0, the Rabi 
frequency is a real parameter, which may be written as 

Qi = Q cos {ki-Ti^ . (104) 

Hence, in the standing-wave laser field the atoms experience different ampli- 
tudes of the driving field. 

In the following, we choose the reference frame such that the atoms are at 
the positions ri = (ri, 0, 0) and r2 = (^^2, 0, 0) along the x-axis, with distance 
ri2 apart. In this case, 

Qi = fie*^^-"^i , ^2 = f^e*^^ "' , (105) 
for the atoms in the running-wave field, and 

rii = cos (kL ■ fi) , ^2 = ^ cos (A;^ ■ r2] , (106) 



for the atoms in the standing-wave field. 

With the above choice of the Rabi frequencies, the Hamiltonian 
takes the form 

Hl = -h{nSt + }i.c.) , (107) 

where = Si exp{ikL ■ fi) + S2 exp{ikL ■ ^2) for the running-wave field, 
and S^ = Si cos^ki ■ ri) + S2 cos{kL ■ ^2) for the standing-wave field. The 
operator S^ corresponds to the symmetric superposition operator defined in 
Eq. (jHH|) . Following the procedure, we developed in Sec. 13. 2. 2^ we find that 
the transformation coefficients u and v are 

for the running-wave field, and 

cos(A;l ■ ri) 



u 



cos^iki ■ n) + cos^iki ■ ^2) 

cos{kL ■ r2) 
cos^iki ■ "^i) + cos^iki ■ r2) 



(109) 



for the standing-wave field. 

Using the transformation coefficients ()108p and ()109p . we find that the 
spontaneously induced coherences Vas and the coherent coupling Ac between 
the symmetric and antisymmetric transitions are 

— # — * 

Vas = -iT sm{kL ■ f^) , Ac = zAi sm{kL ■ ru) , (HO) 
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Figure 6: Time evolution of the total radiation intensity for the running-wave 
driving field with Q = 0.2r, ki \\ n2 and different interatomic separations; 
ri2 = 0.2A (sohd line), = 0.16A (dashed line), = 0.14A (dashed-dotted 
line) . 

for the running-wave field, and 

i as — — i 12 , — "12 , U-LJ-J 

I + cos-^ {kL ■ r 12) 1 + cos^ [kL ■ r 12) 

for the standing-wave field, where, for simplicity, we have chosen the reference 
frame such that ri = and r2 = ri2. 

First, we note that no quantum beats can be obtained for the direction of 
propagation of the laser field perpendicular to the interatomic axis, because 
sm{kL ■ ^12) = 0; however, quantum beats occur for directions of propagation 
different from the perpendicular to fi2. One can see from Eqs. and 
(llllj) that in the case of the running- wave field and Al = 0, the symmetric 
and antisymmetric transitions are correlated only through the spontaneously 
induced coherences Tas- In the case of the standing- wave field, both coupling 
parameters Tas and Ac are different from zero. However, for interatomic 
separations ri2 < A, the parameter Tas is much smaller than A^, indicating 
that in this case the coherent coupling dominates over the spontaneously 
induced coherences. These simple analysis of the parameters Tas and Ac show 
that one should obtain quantum beats in the total radiation intensity of the 
fluorescence field emitted from two identical atoms. Figures IHl and [7| show 
the time evolution of the total radiation intensity, obtained by numerical 
solutions of the equations of motion for the atomic correlation functions. 
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Figure 7: Time evolution of the total radiation intensity for the same param- 
eters as in Fig. IHl but the standing- wave driving field. 

The equations are found from the master equation (j42j) . which in the case 
of the running- or standing-wave driving field leads to a closed set of fifteen 
equations of motion for the atomic correlation functions j5U| I91j. In Fig. IHl 
we present the time-dependent total radiation intensity for the running-wave 
driving field with Q = 0.2r, ki \\ fi2 and different interatomic separations. 
Fig.[7|shows the total radiation intensity for the same parameters as in Fig.lHl 
but the standing- wave driving field. As predicted by Eqs. (jllUj) and (jllip . 
the intensity exhibits quantum beats. The amplitude and frequency of the 
oscillations is dependent on the interatomic interactions and vanishes for 
large interatomic separations as well as for separations very small compared 
with the resonant wavelength. This is easily explained in the framework 
of collective states of a two-atom system. For a weak driving field, the 
population oscillates between the intermediate states |s), \a) and the ground 
state \g). When interatomic separations are large, flu is approximately zero, 
and then the transitions \s) \g) and \a) —>■ \g) have the same frequency. 
Therefore, there are no quantum beats in the emitted field. On the other 
hand, for very small interatomic separations, ■ fi2 ~ 0, and then the 
coupling parameters Tas and vanish, resulting in the disappearance of the 
quantum beats. 
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5 Nonclassical states of light 



The interaction of light with atomic systems can lead to unique phenomena 
such as photon antibunching and squeezing. These effects are examples of a 
nonclassical light field, that is a field for which quantum mechanics is essen- 
tial for its description. Photon antibunching is characteristic of a radiation 
in which the variance of the number of photons is less than the mean number 
of photons, i.e. the photons exhibit sub-Poissonian statistics. Squeezing is 
characteristic of a field with phase-sensitive quantum fluctuations, which in 
one of the two phase components are reduced below the vacuum (shot-noise) 
level. Since photon antibunching and squeezing are distinguishing features 
of light, it is clearly of interest to identify situations in which such fields 
can be generated. Photon antibunching has been predicted theoretically for 
the first time in resonance fluorescence of a two-level atom |1UH llU2j . Since 
then, a number of papers have appeared analyzing various schemes for gen- 
erating photon antibunching offered by nonlinear optics [1031 11U41 11U5| HUB] . 
Squeezing has been extensively studied since the theoretical work by Walls 
and Zoller |lU7j and Mandel |lU8j on reduction of noise and photon statistics 
in resonance fluorescence of a two-level atom. Several experimental groups 
have been successful in producing nonclassical light. Photon antibunching 
has been observed in resonance fluorescence from a dilute atomic beam of 
sodium atoms driven by a coherent laser field |109t IllUL llllj . More re- 
cently, beautiful measurements of photon antibunching have been made on 
trapped atoms and a cavity QED system |113j . On the other hand, 

squeezed light was first observed by Slusher et al. |114j in four-wave mix- 
ing experiments. After that observation squeezed light has been observed in 
many other nonlinear processes, with a recent development being the avail- 
ability of a tunable source of squeezed hght exhibiting a noise reduction of 
~ 70% below the shot-noise level. The experimental observation of photon 
antibunching and squeezing have provided direct evidence of the quantum 
nature of light, and these two phenomena were precursors of much of the 
present work on nonclassical light fields. An extensive literature on various 
aspects of photon antibunching and squeezing now exists and is reviewed in 
several articles jll5t 11161 1117j . 

The objective of this section is to concentrate on collective two-atom 
systems as a potential source for photon antibunching and squeezing. We 
understand collective effects in a broad sense, that for two or more atoms 
all effects that cannot be explained by the properties of individual atoms are 
considered as collective. This definition of collective effects thus includes, 
for example, both the resonance fluorescence from a system of two atoms 
in free space and also collective behaviour of two atoms strongly coupled to 
the same cavity mode in the good cavity limit. Moreover, we emphasize the 
role of the interatomic interactions in the generation of nonclassical light. 
We also relate the nonclassical effects to the degree of entanglement in the 
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system. 



5.1 Photon antibunching 

Photon antibunching is described through the normahzed second-order cor- 
relation function, defined as 



where 

= (E(-) (i?2,t2) (i?2,t2) E^^^ {Ruh)) ,(113) 

is the two-time second-order correlation function of the EM field detected at 
a point Ri at time ti and at a point R2 at time t2, and 

is the first-order correlation function of the field (intensity) detected at a 
point Ri at time ti{i = 1,2). 

The correlation function G^'^\Ri,ti; R2, ^2) is proportional to a joint prob- 
ability of finding one photon around the direction Ri at time ti and another 
photon around the direction R2 at the moment of time ^2- For a coher- 
ent light, the probability of finding a photon around Ri at time ti is in- 
dependent of the probability of finding another photon around R2 at time 
t2, and then G^^\Ri,ti; R2,t2) simply factorizes into G^^\Ri,ti)G^^\R2,t2) 
giving g^'^\Ri,ti; R2,t2) = 1. For a chaotic (thermal) field the second-order 
correlation function for ti = t2 is greater than for ^2 — ^1 = t > giv- 
ing g^'^\Ri,ti; R2,ti) > g^'^\Ri,ti] R2,ti + r). This is a manifestation of 
the tendency of photons to be emitted in correlated pairs, and is called 
photon bunching. Photon antibunching, as the name implies, is the oppo- 
site of bunching, and describes a situation in which fewer photons appear 
close together than further apart. The condition for photon antibunching is 
g^'^\Ri,ti; R2,ti) < g^'^\Ri,ti; R2,ti + r) and implies that the probability 
of detecting two photons at the same time t is smaller than the probability 
of detecting two photons at different times t and t + t. Moreover, the fact 
that there is a small probability of detecting photon pairs with zero time 
separation indicated that the one-time correlation function g^'^\Ri,t; R2,t) 
is smaller than one. This effect is called photon anticorrelation. The nor- 
malized one-time second-order correlation function carries also information 
about photon statistics, which is given by the Mandel's Q parameter defined 
as dUHl 

Q = qT[g^^^[R,,t;R2,t)-l\ , (115) 
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where q is the quantum efficiency of the detector and T is the photon counting 
time. 

We can relate the field correlation functions ()113p and ()114p to the correla- 
tion functions of the atomic operators, which will allow us to apply directly 
the master equation (|42|) to calculate photon antibunching in a collective 
atomic system. The relation between the positive frequency part of the elec- 
tric field operator at a point R = RR, in the far-field zone, and the atomic 
dipole operators 5*^", is given by the well-known expression [17| HU] 



EM{R,t) = Et\R,t) 



ijj- R'x [R'x Ui] / R\ / X 
- y^5r (t - ^) exp . rl) ,(116) 

where uji is the angular frequency of the ith atom located at a point r^, and 
^o"'"'' ^) denotes the positive frequency part of the field in the absence of 
the atoms. 

If we assume that initially the field is in the vacuum state, then the free- 
field part ^o^'' (-R, t) does not contribute to the expectation values of the 
normally ordered operators. Hence, substituting Eq. ()116|) into Eqs. pi3p 
and (|114|) . we obtain 



N 

G^'^ [R,t-R,t + T) = u{R,)u{R2) E (r.r.Tfer 

i,j,k,l=l 



X {St{t)S^{t + r)Sr{t + r)Sr{t)) 

X exp \^ik (^Ri ■ fij + R2 ■ fkij , (117) 

N ^ 

G«(i?,t) = «(i?) E (r.r,)^(5+(t)5r(t)) 

X exp(ikR-fij^ , (118) 



where r = ^2 — ^i, Tj is the damping rate of the ith atom, and u{R) is 
a constant given in Eq. (ffn|) . The second-order correlation function pi7j] 
involves two-time atomic correlation function that can be calculated from 
the master equation (jHUjl or (j^^ and applying the quantum regression theo- 
rem jll8j . From the quantum regression theorem, it is well known that for 
r > the two-time correlation function {S^ (t) {t + r) Sf {t + r) S~ (t)) 
satisfies the same equation of motion as the one-time correlation function 

{s^it)srit)). 

We shall first of all consider the simplest collective system for photon 
antibunching; two identical atoms in the Dicke model. Whilst this model 
is not well satisfied with the present sources of two-atom systems, it does 
enable analytic treatments that allow to understand the role of the collective 
damping in the generation of nonclassical light. 
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For the two-atom Dicke model the master equation (jl^ reduces to 



^ = hn[s+ + S-,p\-h{s+S-p + pS+S--2S-pS+) ,{119) 



where = + S2 and = Sf + are the collective atomic operators 
and Q is the Rabi frequency of the driving field, which in the Dicke model is 
the same for both atoms. For simplicity, the laser frequency ul is taken to 
be exactly equal to the atomic resonant frequency Uq. 

The secular approximation technique has been suggested by Agarwal et 
al. and Kilin |12()j . which greatly simplifies the master equation 
Hassan et al. |121j and Cordes |122| I123j have generalised the method to 
include non-zero detuning of the laser field and the quasistatic dipole-dipole 
potential. The technique is a modification of a collective dressed-atom ap- 
proach developed by Freedhoff |124j and is valid if the Rabi frequency of the 
driving field is much greater than the damping rates of the atoms, i7 ^ F. 
To implement the technique, we transform the collective operators into new 
(dressed) operators 



= ±]-i + R-) + , 



2 

= -U[r^-R-) . (120) 

The operators R are a rotation of the operators S. For a strong driving field, 
the operators R^ vary rapidly with time, approximately as exp{±iQt) , while 
R^ varies slowly in time. By expressing the operators and in terms of 
the operators i?^ and R^, and substituting into the master equation ()119|) . 
we find that certain terms are slowly varying in time while others oscillate 
rapidly. The secular approximation then involves dropping the rapidly oscil- 
lating terms that results in an approximate master equation of the form 

^ = iQ[R',p] - lr{{R'R'p + pR'R' -2R'pR') 



+i (R^R-p + pR+R- - 2R-pR^ 
+j[r-R+P + pR'R+ -2R+pR-)\ . (121) 



4 

The master equation ()12H) enables to obtain equations of motion for the 
expectation value of an arbitrary combination of the transformed operators 
R. In particular, the master equation leads to simple equations of motion 
for the expectation values required to calculate the normalized second-order 
correlation function. The required equations of motion are given by 
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Figure 8: The normalised second-order correlation function g^'^'' (r) as a func- 
tion of r and different fi; f2 = 2.5r (solid line), Q = lOF (dashed line). 



= - (^r + 2.fi) . (122) 

The solution of these decoupled differential equations is straightforward. Per- 
forming the integration and applying the quantum regression theorem |118j , 
we obtain from Eqs. ()122|) and ()112|) the following solution for the normalized 
second-order correlation function jHH I125j 

c/(2)(r) = lim^(^) (i?i,t;^2,t + = 1 + — exp f--rr 



32 ' V 2 

+ - exp [-^Trj cos (2fir) - - exp (-^Fr j cos (^]r) (.123) 

The correlation function g^"^^ (r) is shown in Fig. |H1 as a function of r for 
different Q. For r = 0, the correlation function g^^'' (0) = 0.75, showing 
the photon anticorrelation in the emitted fluorescence field. As r increases, 
the correlation function increases {g^'^^ (r) > g^'^^ (0)), which reflects photon 
antibunching in the emitted field. However, the photon anticorrelation in the 
two-atom fluorescence field is reduced compared to that for a single atom, for 
which g^'^^ (0) = 0. This result indicates that the collective damping reduces 
the photon anticorrelations in the emitted fluorescence field. 

As we have mentioned above, in the Dicke model the dipole-dipole in- 
teraction between the atoms is ignored. This approximation has no justifi- 
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cation, since for small interatomic separations the dipole-dipole parameter 
^12, which varies as {kQri2)~^, is very large and goes to infinity as ri2 goes 
to zero (see Fig. Q). Moreover, the Dicke model does not correspond to the 
experimentally realistic systems in which atoms are separated by distances 
comparable to the resonant wavelength. Ficek et al. [84j and Lawande et 
al. [126J have shown that the dipole-dipole interaction does not considerably 
affect the anticorrelation effect predicted in the Dicke model. Richter |127j 
has shown that the value g^"^"* (0) = 0.75 can in fact be reduced such that 
even the complete photon anticorrelation g^"^"* (0) = can be obtained, if the 
dipole-dipole is included and the laser frequency is detuned from the atomic 
transition frequency. To show this, we calculate the normalized second-order 
correlation function ()112|) for the steady-state fluorescence fleld from two 
identical atoms (A^ = 2), and r = 0. In this case, the correlation func- 
tion (|112|) with Eqs. (|117|) and (jll8|) can be written as 



2f/{l + cos 


kfi2 ■ (-Rl - R2 


)] 


} 


1 + FT cos (kri2 ■ Ri^ 






1 + cos ( 


kfi2 ■ R2^ 





9^'^ (0) 

where U and W are the steady-state atomic correlation functions 

{3^82 Si S2 ) {3^82 + S2 Si ) 



(124) 



U 



{St Si + s+S2y 



w 



{srsi + sts2) 



(125) 



The steady-state correlation functions are easily obtained from the master 
equation (jl^ . We can simplify the solutions assuming that the atoms are 
in equivalent positions in the driving fleld, which can be achieved by propa- 
gating the laser fleld in the direction perpendicular to the interatomic axis. 
In this case we get analytical solutions, otherwise for ki ■ ru 7^ numerical 
methods are more appropriate [201 IHH I128j . With ki ■ ^12 = the master 
equation ()42|) leads to a closed set of nine equations of motion for the atomic 
correlation functions. This set of equations can be solved exactly in the 
steady-state |129j . and the solutions for U and W are 



U 
W 



^4 ^ ^ 4^2 ) ^2 + ^ 4^2 ) h (r + t^)^ + (a 



12 



(12 + 4Ai + 2n^y 



(r^ + 4Ai: 



(12 + 4A2 + 2^]2) 



(126) 



One can see from Eqs. (|124|) and (|126|) that there are two different processes 
which can lead to the total anticorrelation, (7*^2) (^q^ = q ^j^g g^-g^ one involves 
an observation of the fluorescence fleld with two detectors located at different 
points. If the correlation function is measured using two detectors, Ri R2, 
and then we obtain g^"^^ (0) = whenever the positions of the detectors are 
such that 



{l + cos [kri2 ■ [R1-R2)]} = 



(127) 
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which happens when 



kfu ■ {Ri - R2) = {2n + 1) 



71 



n = 0,±1,±2. . . 



(128) 



In other words, two photons can never be simultaneously detected at two 
points separated by an odd number of A/2ri2, despite the fact that one 
photon can be detected anywhere. This complete anticorrelation effect is 
due to spatial interference between different photons and reflects the fact 
that one photon must have come from one source and one from the other, 
but we cannot tell which came from which. 

It should be emphasised that this effect is independent of the interatomic 
interactions and the Rabi frequency of the driving field. The vanishing of 
g^"^^ (0) for two photons at widely separated points ^1 and R2 is an example 
of quantum-mechanical nonlocality, that the outcome of a detection mea- 
surement at -Ri appears to be influenced by where we have chosen to locate 
the i?2 detector. At certain positions R2 we can never detect a photon at 
Ri when there is a photon detected at R2, whereas at other position R2 it 
is possible. The photon correlation argument shows clearly that quantum 
theory does not in general describe an objective physical reality independent 
of observation. 

The second process involves the shift of the collective atomic states due 
to the dipole-dipole interaction that can lead to ^(2) (0) = even if the 
correlation function is measured with a single detector (^1 = ^2) or two 
detectors in configurations different from that given by Eq. (jl28j) . For a 
weak driving field (f2 <C F) and large detunings such that Al = Q12 ^ F, 
the correlation function (|112j) with ^1 = R2 simplifies to 



Thus, a pronounced photon anticorrelation, gf*^^^ (0) ~ 0, can be obtained for 
large detunings such that = Q12, i.e., when the dipole-dipole interaction 
shift of the collective states and the detuning cancel out mutually. The 
correlation function g^'^^ (0) of the steady-state fluorescence field is illustrated 
graphically in Fig. IHlas a function of A^, for the single detector configuration 
with Ri = R2 = R, and different ru- The graphs show that g^'^^ (0) strongly 
depends on A^, and the total photon anticorrelation can be obtained for 
A^, = fli2- Referring to Fig. the condition A^, = flu corresponds to the 
laser frequency tuned to the resonance with the \g) — > \s) transition. Since 
the other levels are far from the resonance, the two-atom system behaves like 
a single two-level system with the ground state \g) and the excited state \s). 
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(2) (0) 



(F + F12)' 
4Ai 



(129) 
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Figure 9: The normalised second-order correlation function g^'^'> (0) as a func- 
tion of Al for Ri = R2 = R,ri2 -L R,p, -L fu, fl = 0.5T and different ri2; 
ri2 = lOA (solid line), — 0.15A (dashed line), — 0.08A (dashed-dotted 
line) . 

5.2 Squeezing 

To understand squeezed light, recall that the electric field amplitude E (f) 
may be expressed by positive- and negative-frequency parts 

E{r)^ ^(+) (r) + (r) , (130) 

where 

Ei+) (r-) = (^(-) {r})' = -i ^ {hco,/2eoVf' ^^s^^'^^ , (131) 



ks 



and ujk = c 



is the angular frequency of the mode k. 
We introduce two Hermitian combinations (quadrature components) of 
the field components that are 7r/2 out of phase as 



-ie 



Ee = (i?) e^^ + (i?) e- 

Ee-./2 = -^ {e^^^ (r) e^' - E^-^ {r) e-'^') , (132) 

where 

e^ut-k-R, (133) 



45 



and uj is the angular frequency of the quadrature components. 

The quadrature components do not commute, satisfying the commutation 
relation 

where C is a positive number 

C = Y.\rii^kl'^eoV\ . (135) 

fcs 

Hence the two quadrature components cannot be simultaneously precisely 
measured, and from the Heisenberg uncertainty principle, we find that the 
variances ^A^|^ and {j\.El_^i^ satisfy the inequality 

(a4^) (A^^^/,) > , (136) 

where the equality holds for a minimum uncertainty state of the field. 

The variances l^lS.E'^ and {jS.El_^i^ depend on the state of the field 
and can be larger or smaller than C . A chaotic state of the field leads to the 
variances in both components larger than C: 

{/\El)>C and (AE^^/2)>C. (137) 

If the field is in a coherent or vacuum state 

(A^^) = {/\El_^i,) = C , (138) 

which is an example of a minimum uncertainty state. 

A squeezed state of the field is defined to be one in which the variance in 
one of the two quadrature components is less than that for the vacuum field 

(AE|)<C or {/\El_^,,)<C . (139) 

The variances can be expressed as 

(A^^) = C+{:/\El-) , 

(Ai|_,/2> = Ai^^/^:) , (140) 

where the colon stands for normal ordering of the operators. 

As the squeezed state has been defined by the requirement that either 
(^AEq'^ or (^AEl_^i^ be below the vacuum level C, it follows immediately 
from Eq. ()140j) that either 

(:A42:)<0 or (: AE^^/^ :) < (141) 
for the field in a squeezed state. 



2iC 



(134) 
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We now determine the relation between variances in the field and the 
atomic dipole operators. Using Eq. ()116p . which relates the field operators 
to the atomic dipole operators, we obtain 



u R 



(142) 



where a = 9,9 — 11/2, Sa and 5*3 are real (phase) operators defined as 

Se = I (5+ + S,) , = ^ {S^ - S,) , (143) 



and 



with 



1 



Sa^ , Sfj 



N 



i=l 



±i[kR- f,- - 9 



(144) 



(145) 



We first consider quantum fiuctuations in the fiuorescence field emitted 
by two identical atoms in the Dicke model. To simplify the calculations we 
will treat only the case of zero detuning, A/^ = 0. Assuming that initially 
(t = 0) the atoms were in their ground states, we find from Eqs. ()122|) and 
the following expressions for the time-dependent variances |13Uj 



i-iexp (^-^n)cos(2fit) 



Fg=,{t) ^ {: AE',^, :) / {2u{R) 



1 / 3 

- exp f — -Ft ) cos (fit) 



(146) 



and 



'e=Tv/2 



it) = 



■■)/{mr) 



-exp 



--Tt] cos (nt) 



(147) 



In writing Eqs. (jl46|) and (jl47|) . we have assumed that the angular frequency 
of the quadrature components is equal to the laser frequency, u = ul, and 
we have normalised the variances such that F{t) determines fiuctuations per 
atom. It is easily to show that the variance Fg=T^/2 (t) is positive for all 
times t, and squeezing {Fq < 0) can be observed in the variance Fg=Q{t). 
The time dependence of the variance Fg=Q (t) is shown in Fig. El for two 
different values of the Rabi frequency. It is seen that squeezing appears in the 
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Figure 10: The variance Fg=Q (t) as a function of time for different f2; f2 = 
lOOr (solid line), n = 200r (dashed line). 

transient regime of resonance fluorescence and its maximum value (minimum 
of F) moves towards shorter times as Q increases. The optimum squeezing 
reaches a value of —1/16 at a very short time. This value is equal to the 
maximum possible squeezing in a single two-level atom |1U7| I13H I132t 1133] . 
Thus, the collective damping does not affect squeezing in the two-atom Dicke 
model. This is in contrast to the photon anticorrelation effect which is greatly 
reduced by the collective damping. 

Figure ^1 shows that in the two-atom Dicke model there is no squeezing 
in the steady-state resonance fluorescence when the atoms are excited by 
a strong laser field. Ficek et al. |129j and Richter |134j have shown that 
similarly as in the case of photon anticorrelations, a large squeezing can be 
obtained in the steady-state resonance fluorescence from a strongly driven 
two-atom system, if the dipole-dipole interaction is included and the laser 
frequency is detuned from the atomic transition frequency. This is shown 
in Fig. El where we plot -Fe=0; calculated from Eq. ()142|) and the master 
equation (j^ . for the steady-state resonance fluorescence from two identical 
atoms, with fi2 _L -R, f2 = O.SF, ■ "^12 = and different ri2. It is evident 
from Fig. ^2 that a large squeezing can be obtained for a finite A/, and its 
maximum shifts towards larger A^;, as the interatomic separation decreases. 
Similar as in the case of photon anticorrelations, the maximum squeezing 
appears at A^ = f2i2, and can again be attributed to the shift of the collective 
energy states due to the dipole-dipole interaction. 
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Figure 11: The steady-state variance Fq=q as a function of for fi2 -L R, 
fi _L fi2, f2 = O.Sr, ■ ri2 = and different ru; ri2 = lOA (solid line), 
ri2 = 0.15A (dashed line), ri2 = 0.08A (dashed-dotted line). 

The variance -fe=05 shown in Fig. lll[ exhibits not only the large squeezing 
at finite detuning A^, but also a small squeezing near = 0. In contrast 
to the squeezing at finite A^, which has a clear physical interpretation, the 
source of squeezing at Al is not easy to understand. To find the source of 
squeezing at A^ = 0, we simplify the calculations assuming that the angular 
frequency of the quadrature components u = and the fluorescence field 
is observed in the direction perpendicular to the interatomic axis, R -L fi2. 
In this case, the variance AE'^ written in terms of the density matrix 
elements of the collective system, is given by |135j 

F^ = {:AEl:)/(2u{R)) = 1 {2pee + 2p,, + pe.e^^'^ + p.ee"'^" 

- '{Pes + Psg) + {Pse + Pgs) e"^"]'} (148) 

This equation shows that the variance depends on phase a not only through 
the one-photon coherences Pes and p^^, but also through the two-photon 
coherences p^g and Pge- This dependence suggests that there are two different 
processes that can lead to squeezing in the two-atom system. The one- 
photon coherences cause squeezing near one-photon resonances |e) — * \s) 
and \s) \g), whereas the two-photon coherences cause squeezing near the 
two-photon resonance \g) \e). 
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To show this, we calculate the steady-state populations and coherences 
from the master equation (jl^ . We use the set of the collective states (jHSI) 
as an appropriate representation for the density operator 



P 



(149) 



where pij are the density matrix elements in the basis of the collective states. 

After transforming to the collective state basis, the master equation ()42p 
leads to a closed system of fifteen equations of motion for the density matrix 
elements jHl]. However, for a specifically chosen geometry for the driving 
field, namely that the field is propagated perpendicularly to the atomic axis 
i^L ■ = 0), the system of equations decouples into nine equations for 
symmetric and six equations for antisymmetric combinations of the density 
matrix elements jHH ESI HOI HH 11^71 1128] . In this case, we can solve the 
system analytically, and find that the steady-state values of the populations 
and coherences are I127j 
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Psg 



Peg 
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in^ (r + 2iAL) /z , 

n [m (n + 2iAL 



^p2 ^ 4^2 ) ^ ^4 



+ (r^ + 4A| 



(r + 2iAi) 



^ (r + ri2) + ^ (Ai - ^^12) 



/z 
/z, 



where 



z = 4fi^ + (r2 + 4A2) \ 2n^ + 



- (r + ri2)' + {Al-q 



12 



(150) 



(151) 



and n = n/V2. 

Near the one-photon resonance \s) \g) the detuning Al = ^u, and 
assuming that fli2 ^ fl, T, the coherences reduce to 



Pe 



12 



Psg 



-2 (r + r 



12 
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i{T + T 



12 



12 



(152) 



It is clear from Eq. p52|l that near the \g) \s) resonance the coherence 
Psg is large, whereas the two-photon coherence is of order of ili2 and thus is 
negligible for large Qu- 

Near two-photon resonance, A^ ^ 0, and it follows from Eq. ()150|) that 
in the limit of Q12 ^ Q,T the coherences reduce to 



pe 



psg 



n 



12 



peg 



(153) 
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Figure 12: The steady-state variance Fa as a function of A/^ for ri2 = 0.05A, 
i7 = 3r, ^ _L ri2, /i J- fi2 and different phases a: a = tx/2 (dashed hne), 
a = 3-71 / 4 (sohd hne). 



In this regime, the coherences psg and p^g are of order of magnitude , 
but Peg dominates over the one-photon coherence Psg when the driving field 
is strong |136j . 

The steady-state variance Fa, calculated from Eqs. ()148p and pSOj) . is 
plotted in Fig. [T21 as a function of for ri2 = 0.05 A, fl = 3T, R ± ri2 
and different phases a. The variance shows a strong dependence on a near 
the one- and two-photon resonances. Moreover, a large squeezing is found 
at these resonances. It is also seen that near the two-photon resonance a 
change by 7r/4 of the phase a changes a dispersion-like structure of Fa into 
an absorption-hke type. According to Eqs. p48|) and p50|) . the variance Fa 
for Q12 ^ S> r, can be written as 



^12 



Al r 

cos 2a + ——: — TTTT sin 2a 



(r2 + 4Ai) (r2 + 4A2 



L) 



(154) 



where we retained only those terms which contribute near the two-photon 
resonance. Equation ()154|) predicts a dispersion-like structure for a = or 
7^/2, and an absorption-like structure for a = n/A. Moreover, we see that 
the presence of the dipole-dipole interaction is essential to obtain squeezing 
near the two-photon resonance. The emergence of an additional dipole-dipole 
interaction induced squeezing is a clear indication of a totally different process 
which can appear in a two-atom system. The dipole-dipole interaction shifts 
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the collective states that induces two-photon transitions responsible for the 
origin of the two-photon coherence. 



6 Quantum interference of optical fields 

In the classical theory of optical interference the EM field is represented by 
complex vectorial amplitudes E{R, t) and E*{R, t), and the first- and second- 
order correlation functions are defined in a similar way as the correlation 
functions (ITT!?|l and (fTTl) of the field operators E^+\R,t) and E^-\R,t). 
This could suggest that the only difference between the classical and quan- 
tum correlation functions is the classical amplitudes E*{R,t) and E{R,t) 
are replaced by the field operators E^~''{R,t) and E^^\R,t). This is true 
as long as the first-order correlation functions (coherences) are considered, 
where the interference effects do not distinguish between the quantum and 
classical theories of the EM field |ll-{7j . However, there are significant dif- 
ferences between the classical and quantum descriptions of the field in the 
properties of the second-order correlation function jHEl 1138] . 



6.1 First-order interference 

The simplest system in which the first-order interference can be demonstrated 
is the Young's double slit experiment in which two light beams of amplitudes 
Ei{fi,ti) and E2{r2,t2), produced at two slits located at ri and r2, respec- 
tively, incident on a detector located at a point R far away from the slits. 
The resulting average intensity of the two fields measured by the detector 
can be written as 



{I{R,t)) = a{{h{Ri,t-h)) + {l2{R2,t~t2)) 

+2Re{E*{Rut-U)E2{R2,t-t2))} , (155) 

where a is a constant that depends on the geometry and the size of the slits, 
and Ii{Ru t-U) = {E*{Ri, t - U)Ei{Ri, t - U)). 

If the observation point R lies in the far field zone of the radiation emitted 
by the slits, the fields at the observation point can be approximated by plane 
waves for which we can write 

Ei{Ri,t-ti) ^ Ei{R,t)exp[-i(uit-kiR-fi + (j)i^^ , (156) 

where ki = uji/c, uoi is the angular frequency of the iih field and 0j is its 
initial phase. 

For perfectly correlated fields with equal amplitudes and frequencies, and 
fixed the phase difference 0i — 02, the average intensity detected at the 
point R is given by 

{I{R,t)) = 2<y{h){l + cos kR-r 12) , (157) 
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where Iq = Ii = I2. 

Equation (jlSTj) shows that the average intensity depends on the position 
R of the detector, and small changes in the position R of the detector lead 
to minima and maxima in the detected intensity. The usual measure of 
the minima and maxima of the intensity, called the interference fringes or 
interference pattern, is a visibility defined as 

•> 1 I max Imin / 1 r o \ 

V = , (158) 

-'■max ~ -'■mm 

where Imax corresponds to cos{kR ■ fi2) = 1, whereas Imin corresponds to 
cos{kR-fi2) = —1 of the field intensity ()157|) . The visibility of the interference 
fringes corresponds to the degree of coherence between two fields. Hence, 
two classical fields of equal amplitudes and frequencies, and fixed the phase 
difference produce maximum possible interference pattern with the maximum 
visibility of 100%. 

For a quantum field, the electric field components can be expressed in 
terms of plane waves as 

(f, t) = (E(-) if, t))' = E (1^) ' %,a,,e<^^---*) , (159) 

ks ^ ' 

where V is the volume occupied by the field, a^^ is the annihilation operator 
for the fcth mode of the field of the polarization e^^ and 0;^ is the angular 
frequency of the mode. 

It is easily to show, that in the case of interference of quantum fields, 
the average intensity detected at the point R has the same form as for the 
classical fields, Eq. ()157|) . with (/q) given by 

(^0) = E^(-J, (160) 

ks 

J, — * 

where (^^3) = {cb^^J^ks) average number of photons in the mode k. 

Thus, interference effects involving the first-order coherences cannot dis- 
tinguish between the quantum and classical theories of the EM field. 

6.2 Second-order interference 

The second-order correlation function has completely different coherence 
properties than the first-order correlation function. An interference pattern 
can be observed in the second-order correlation function even if the fields are 
produced by two independent sources for which the phase difference 0i — 02 
is completely random |139| I140j . In this case the second-order correlation 
function, observed at two points Ri and i?2, is given by |138j 

G^'\R,,t,;R2,t2) = {Il{h)) + {ll{t2))+2{h{t,)l2{t2)) 

+ 2(/i(ti)/2(t2))cos[A;fi2- (i?i-i?2)l . (161) 
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Clearly, the second-order correlation function of two independent fields ex- 
hibits a cosine modulation with the separation Ri — R2 of the two detectors. 
This is an interference although it involves a correlation function that is of 
the second order in the intensity. Similar to the first-order correlation func- 
tion, the sharpness of the fringes depends on the relative intensities of the 
fields. For classical fields of equal intensities, Ii = I2 = lo, the correlation 
function 1)1611) reduces to 

G(2)(i?i,t;i?2,t) =4(j2)|l + icos[A;fi2- (i?i-/22)]} . (162) 

In analogy to the visibility in the first-order correlation function, we can 
define the visibility of the interference pattern of the intensity correlations as 

y(2) ^ ^ma. ^rnin ^^gg^ 



" mm 



and find from Eq. ()162|) that in the case of classical fields an interference 
pattern can be observed with the maximum possible visibility of V*-^-* = 1/2. 
Thus, two independent fields of random and uncorrelated phases can exhibit 
an interference pattern in the intensity correlation with a maximum visibility 
of 50%. 

As an example of second-order interference with quantum fields, consider 
the simple case of two single-mode fields of equal frequencies and polariza- 
tions. Suppose that there are initially n photons in the field Ei and m 
photons in the field i?2, and the state vectors of the fields are the Fock states 
= \n) and \ip2) = |^)- The initial state of the two fields is the direct 
product of the single-field states, \ip) = \n)\m). Inserting Eq. ()159|) into 
Eq. ()113|) and taking the expectation value with respect to the initial state 
of the fields, we find 

{n{n-l)+m{m-l) 

+2nm [1 + cos kfu ■ - R2)] } . (164) 

We note that the first two terms on the right-hand side of Eq. p64|l vanish 
when the number of photons in each field is smaller than 2, i.e. n < 2 and 
m < 2. In this limit the correlation function (11641) reduces to 



(165) 



h; R2, ts) = 2 [ ^ j [1 + cos kfu ■ {Ri - R 



Thus, perfect interference pattern with the visibility V*-^-* = 1 can be observed 
in the second-order correlation function of two quantum fields each containing 
only one photon. According to Eq. p62j) . the classical theory predicts only 
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a visibility of V*-^-* = 0.5. For n, m ^ 1, the first two terms on the right-hand 
side of Eq. ()164|) are different from zero (m(m — 1) ^ n{n — 1) ^ n"^), and 
then the quantum correlation function ()164|) reduces to that of the classical 
field. 

The visibility of the interference pattern of the intensity correlations pro- 
vides a means of testing for quantum correlations between two light fields. 
Mandel et al. |141l 11421 1143j have measured the visibility in the interference 
of signal and idler modes simultaneously generated in the process of degen- 
erate parametric down conversion, and observed a visibility of about 75%, 
that is a clear violation of the upper bound of 50% allowed by classical cor- 
relations. Richter |144j have extended the analysis of the visibility into the 
third-order correlation function, and have also found significant differences in 
the visibility of the interference pattern of the classical and quantum fields. 



6.3 Quantum interference in two- atom systems 

In the Young's interference experiment the slits can be replaced by two atoms 
and interference effects can be observed between coherent or incoherent fields 
emitted from the atoms. The advantage of using atoms instead of slits is that 
a given time each atom cannot emit more than one photon. Therefore, the 
atoms can be regarded as sources of single photon fields. 
Using Eq. ()116p . we can write the visibility as 

V = , (166) 

(Si Si + 5*2 S2 ) 

which shows that the interference effects can be studied in terms of the atomic 
correlation functions. 

There have been several theoretical studies of the fringe visibility in 
the fluorescence field emitted by two coupled atoms |145j . and the Young's 
interference-type pattern has recently been observed experimentally in the 
resonance fluorescence of two trapped ions The experimental results 

have been explained theoretically by Wong et al. |146j . and can be under- 
stood by treating the ions as independent radiators which are synchronized 
by the constant phase of the driving field. It has been shown that for a weak 
driving field, the fluorescence field is predominantly composed of an elastic 
component and therefore the ions behave as point sources of coherent light 
producing an interference pattern. Under strong excitation the fluorescence 
field is mostly composed of the incoherent part and consequently there is no 
interference pattern. To show this, we consider a two-atom system driven by 
a coherent laser field propagating in the direction perpendicular to the in- 
teratomic axis. In this case, we can use the master equation (j^^ and obtain 
the analytical formula for the fringe visibility of the steady-state fluorescence 
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field as (HII 



V 



(r^ + 4Ai 



(P + 4Ai) 



(167) 



It is seen that in this specific case, the visibility is positive for all param- 
eter values and is independent of the interatomic interactions. For a weak 
driving field, VL <^ T,Al, the fringe visibility \V\ ~ 1, whereas \V\ ~ 
for VL S> r. A/,, showing that and there is no interference pattern when the 
atoms are driven by a strong field. For moderate Rabi frequencies, ^ F, 
the visibility may be improved by detuning the laser field from the atomic 
resonance. Kochan et al. |148j have shown that the interference pattern of 
the strongly driven atoms can also be improved by placing the atoms inside 
an optical cavity. The coupling of the atoms to the cavity mode induces 
atomic correlations which improves the fringe visibility. 

Here, we derive general criteria for the first- and second-order interfer- 
ence in the fluorescence field emitted from two two-level atoms. Using these 
criteria, we can easily predict conditions for quantum interference in the two 
atom system. In this approach, we apply the collective states of a two-atom 
system, and write the atomic correlation functions in terms of the density 
matrix elements of the collective system as 



^2 ^2 



Ps 



Paa + 2p 

ee ) 



Si 5*2 

Q+ Q+ Q— Q- 
Di O2 ^1 '-'2 



- {Pss - Paa + pas Psa) 
Pee ! 



(168) 



where pu^i = a,s,e) are the populations of the collective states and Psa, Pas 
are coherences. 

From the relations ()168|) . we find that in terms of the density matrix 
elements the first-order correlation function can be written as 



(^R^ t) = Tu{R) {2pee (t) + pss (t) (l + COS kR 
+ Paa (t) (1 - COS kR ■ fis) 

+ i {Psa (t) - Pas {t)) sin kR ■ fi2} , 
and the second-order correlation function takes the form 
G^'^URi,t-R2,t] = AT\{Ri)u{R2) 



■ ri2 



(169) 



X pee (t) 



1 + COS 



k(^Ri 



R, 



ri2 



;i7o) 



It is evident from Eq. ()169p that first-order correlation function can exhibit an 
interference pattern only if p^^ 7^ Paa and/or lm{psa) 7^ 0. This happens when 
{^i\{92\p\^2)\gi) and (5'i|(e2|p|(y'2)|ei) are different from zero, i.e. when there 
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are nonzero coherences between the atoms. Schon and Beige |H21 have arrived 
to the same conclusion using the quantum jump method. On the other hand, 
the second-order correlation function is independent of the populations of 
the entangled states Pss,Paa and the coherences, and exhibit an interference 
pattern when Pee(t) 7^ 0. 

We now examine some specific processes in which one can create unequal 
populations of the \s) and \a) states. Dung and Ujihara |149j have shown 
that spontaneous emission from two identical atoms, with initially only one 
atom excited, can exhibit an interference pattern. Their results can be easily 
interpreted in terms of the populations Pss(t) and Paaif)- If initially only one 
atom was excited; Pee(O) = and Pss(O) = Paa(O) = Psa(O) = Pas(O) = \. 
Using the master equation with Vti = = 0, we find that the time 
evolution of the populations Pssif) and Paaif) is given by 



Since the populations decay with different rates, the symmetric state decays 
with an enhanced rate F + ri2, while the antisymmetric state decays with a 
reduced rate V — ri2, the populations Paaif) is larger than Pssif) for all t > 0. 
Hence, an interference pattern can be observed for t > 0. This effect arises 
from the presence of the interatomic interactions (ri2 7^ 0). Thus, for two 
independent atoms the populations decay with the same rate resulting in the 
disappearance of the interference pattern. 

When the atoms are driven by a coherent laser field, an interference pat- 
tern can be observed even in the absence of the interatomic interactions. To 
show this, we consider the steady-state solutions p50|) for the populations 
of the collective atomic states. It is evident from Eq. that pss > paa 

even in the absence of the interatomic interactions (ri2 = ^^12 = 0). Hence, 
an interference pattern can be observed even for two independent atoms. In 
this case the interference pattern results from the coherent synchronization 
of the oscillations of the atoms by the constant coherent phase of the driving 
laser field. 

We have shown that the first-order coherence is sensitive to the inter- 
atomic interactions and the excitation field. In contrast, the second-order 
correlation function can exhibit an interference pattern independent of the 
interatomic interactions and the excitation process |15mil5l]ll52j . According 
to Eq. ()170|) . to observe an interference pattern in the second-order correla- 
tion function, it is enough to produce a non-zero population in the state 
|e). The interference results from the detection process that a detector does 
not distinguish between two simultaneously detected photons. As an exam- 
ple, consider spontaneous emission from two identical and also nonidentical 
atoms with initially both atoms excited. 

For two identical atoms, we can easily find from Eqs. and (|17(jp . 
and the quantum regression theorem |118j . that the two-time second-order 



Pss if) 

Paa (t) 



p,, (0) exp [- (r + ri2) t] , 
p„, (0) exp [- (r - ri2) t] . 



(171) 
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correlation function is given by 



G(') + = -r\{Ri)u{R2)exp[-r{2t + T)] 



X 



I 1 + cos (j^Ri ■ fi2^ COS (j£R2 ■ ri2^ cosh (ri2r 



cos 



[kR 



1 ■ ri2 



cos 



{kR 



2 ■ ri2 



sinh (ri2r) 



+ sin ■ ris) sin (A;i?2 ■ ru) cos (2ni2r)} . (172) 

The above equation shows that the two-time second-order correlation func- 
tion exhibits a sinusoidal modulation in space and time. This modulation 
can be interpreted both in terms of interference fringes and quantum beats. 
The frequency of quantum beats is 2Qi2 and the amplitude of these beats 
depends on the direction of observation in respect to the interatomic axis. 
The quantum beats vanish for directions 6i = 90° or 62 = 90°, where ^1(^2) 
is the angle between ri2 and Ri{R2), and the amplitude of the beats has 
its maximum for two photons detected in the direction 61 = 62 = 0°. This 
directional effect is connected with the fact that the antisymmetric state \a) 
does not radiate in the direction perpendicular to the interatomic axis. We 
will discuss this directional effect in more details in Sec. 18.11 For independent 
atoms, ri2 = 0, Q12 = 0, and then the correlation function (jl72|l reduces to 



G^^^ [Ri,t;R2,t + T) = -T\{Ri)u{R2)exp[-T{2t + T)] 



X 



1 + cos A; (^Ri - _R2) ■ n2 



(173) 



which shows that the time modulation vanishes. This implies that quantum 
beats are absent in spontaneous emission from two independent atoms, but 
the spatial modulation is still present. 

The situation is different for two nonidentical atoms. In this case, the 
two-time second-order correlation function exhibits quantum beats even if the 
atoms are independent. For ri2 = and ^12 = 0, the master equation (ji^ 
leads to the following correlation function 



G'(') (^i,i;^2,t + r) = -TMRiHR2)exp[~T{2t + T)] 



X 



+ cos 



jcosh^ (r2 - Ti)t 

k - R2) ■ fi2 - 2 At 



} • (174) 



Thus, for independent nonidentical atoms, the correlation function shows a 
sinusoidal modulation both in space and time. We note that the modula- 
tion term in Eq. ()174p is the same as that obtained by Mandel |153j . who 
considered the second-order correlation function for two beams emitted by 
independent lasers. 
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7 Selective excitation of the collective atomic 
states 



In the previous section, we have shown that nonclassical effects in coher- 
ently driven two-atom systems reffect the preparation of the system in a 
superposition of two collective states. In particular, for the total photon 
anticorrelation and maximum squeezing, the two-atom system is in a super- 
position of the ground and the entangled symmetric states. The other states 
are not populated. We now consider excitation processes which can lead to 
a preparation of the two-atom system in only one of the collective states. In 
particular, we will focus on processes which can prepare the two-atom system 
in the entangled symmetric state |s). Our main interest, however, is in the 
preparation of the system in the maximally entangled antisymmetric state 
I a) which, under the condition = y/TiT^, is a decoherence-free state. The 
central idea is to choose the distance between the atoms such that the re- 
sulting level shift is large enough to consider the possible transitions between 
the collective states separately. This will allow to make a selective excitation 
of the symmetric and antisymmetric states and therefore to create controlled 
entanglement between the atoms. 

7.1 Preparation of the symmetric state by a pulse laser 

Beige et al. |154j have shown that a system of two identical two-level atoms 
may be prepared in the symmetric state \s) by a short laser pulse. The 
conditions for a selective excitation of the collective atomic states can be 
analyzed from the interaction Hamiltonian of the laser field with the two- 
atom system. We make the unitary transformation 

Hl = e'^'^^'^HLe-'^^'/^ , (175) 

where 

Ha = h {^L{\e){e\-\g){g\) + {^L + ^12) \s){s\ 

+ {^L-n^2)\a){a\} , (176) 

and find that in the case of identical atoms, Fi = r2 and A = 0, the trans- 
formed interaction Hamiltonian is given by 

Hl = -^{(l]l+f^2)(^eV^''^''''"^* + '^sy^''""''"^*) 

+ (l]2-f^i)(5+e'(^^+^^^)* + 5+e*(^^-^^^)*)+H.c.} . (177) 

The Hamiltonian p77|l represents the interaction of the laser field with the 
collective two-atom system, and in the transformed form contains terms 
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oscillating at frequencies (A^, ± ^12), which correspond to the two sepa- 
rate groups of transitions between the collective atomic states at frequencies 
ul = Uq + Q12 and = ujq — Qi2- The + Q12 frequencies are separated 
from Al — ^12 frequencies by 2Qi2, and hence the two groups of the tran- 
sitions evolve separately when Q12 ^ F. Depending on the frequency, the 
laser can be selectively tuned to one of the two groups of the transitions. 
When ul = ujq + Qu (A^ — Q12 = 0) the laser is tuned to exact resonance 
with the |e) — \a) and \g) — \s) transitions, and then the terms, appearing 
in the Hamiltonian (jl77p . and corresponding to these transitions have no 
explicit time dependence. In contrast, the \g) — \a) and |e) — |s) transitions 
are off-resonant and the terms corresponding to these transitions have an 
explicit time dependence exp(±2if2i2t). If ^^12 ^ T, the off-resonant terms 
rapidly oscillate with the frequency 2Qi2, and then we can make a secular 
approximation in which we neglect all those rapidly oscillating terms. The 
interaction Hamiltonian can then be written in the simplified form 



Hr 



2V2 



(r]i + fi2)5+ + (^^2-^^i)5+ +H.C. . (178) 



It is seen that the laser field couples to the transitions with significantly 
different Rabi frequencies. The coupling strength of the laser to the \g) — \ s) 
transition is proportional to the sum of the Rabi frequencies Q1+Q2, whereas 
the coupling strength of the laser to the \a) — \e) transition is proportional 
to the difference of the Rabi frequencies Qi — Q2- According to Eq. 
the Rabi frequencies Qi and ^^2 of two identical atoms differ only by the 
phase factor exp^iki ■ r^)- Thus, in order to selectively excite the \g) — \s) 
transition, the driving laser field should be in phase with both atoms, i.e. 
Qi = 0,2- This can be achieved by choosing the propagation vector ki of 
the laser orthogonal to the line joining the atoms. Under this condition we 
can make a further simplification and truncate the state vector of the system 
into two states \g) and |s). In this two-state approximation we find from the 
Schrodinger equation the time evolution of the population Ps(t) of the state 
\s) as 

P,(t) = sin^ i^^t] , (179) 



where Q = Qi = ^2- 

The population oscillates with the Rabi frequency of the |s) — \g) tran- 
sition and at certain times Ps{t) = 1 indicating that all the population is in 
the symmetric state. This happens at times 

T„ = (2n+1)-|^, n = 0,l,.... (180) 

Hence, the system can be prepared in the state \s) by simply applying a laser 
pulse, for example, with the duration Tq, that is a standard vr pulse. 
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The two-state approximation is of course an idealization, and a possibility 
that all the transitions can be driven by the laser imposes significant limits 
on the Rabi frequency and the duration of the pulse. Namely, the Rabi 
frequency cannot be too strong in order to avoid the coupling of the laser to 
the |s) — |e) transition, which could lead to a slight pumping of the population 
to the state |e). On the other hand, the Rabi frequency cannot be too small 
as for a small Q the duration of the pulse, required for the complete transfer 
of the population into the state \s), becomes longer and then spontaneous 
emission can occur during the excitation process. Therefore, the transfer 
of the population to the state \s) cannot be made arbitrarily fast and, in 
addition, requires a careful estimation of the optimal Rabi frequency, which 
could be difficult to achieve in a real experimental situation. 

7.2 Preparation of the antisymmetric state 
7.2.1 Pulse laser 

If we choose the laser frequency such that A^, + Q12 = 0, the laser field is 
then resonant to the \a) — \g) and |e) — \s) transitions and, after the secular 
approximation, the Hamiltonian ()177p reduces to 



Hr 



h 



2V2 



{n2-n^)s+ + {n^ + n2)s+ + R.c.] . (isi) 



Clearly, for Qi = —Q2 the laser couples only to the \a) — \g) transition. 
Thus, in order to selectively excite the \g) — \a) transition, the atoms should 
experience opposite phases of the laser field. This can be achieved by choosing 
the propagation vector ki of the laser along the interatomic axis, and the 
atomic separations such that 

ki ■ fi2 = (2n + 1) vr, n = 0, 1, 2, . . . , (182) 

which corresponds to a situation that the atoms are separated by a distance 
ri2 = {2n + l)A/2. 

The smallest distance at which the atoms could experience opposite phases 
corresponds to ri2 = A/2. However, at this particular separation the dipole- 
dipole interaction parameter Qu is small, see Fig. ^ and then all of the 
transitions between the collective states occur at approximately the same 
frequency. In this case the secular approximation is not valid and we cannot 
separate the transitions at A^, + ^12 from the transitions at A^ — Q12. 

One possible solution to the problem of the selective excitation with op- 
posite phases is to use a standing laser field instead of the running-wave 
field. If the laser amplitudes differ by the sign, i.e. E^^ = —E^^ = Eq, and 
ki^ ■ ri = — /cl2 ■ ^2, the Rabi frequencies experienced by the atoms are 

2i ^ / 1 ^ 
^^1 = -^/Lti ■ ^0 sin {-kL ■ fi2 
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where = k^^ = k^^ and we have chosen the reference frame such that 
ri = |ri2 and r2 = ~|^i2- It follows from Eq. ()183p that the Rabi frequen- 
cies oscillate with opposite phases independent of the separation between 
the atoms. However, the magnitude of the Rabi frequencies decreases with 
decreasing r^. 

7.2.2 Indirect driving through the symmetric state 

We now turn to the situation of non-identical atoms and consider different 
possible processes of the population transfer to the antisymmetric state which 
could be present even if the antisymmetric state does not decay to the ground 
level. This can happen when = ^/TiT^; i-e. when the separation between 
the atoms is negligible small. Under this condition the antisymmetric state 
is also decoupled from the driving field. According to Eq. ()95|) . the antisym- 
metric state can still be coupled, through the coherent interaction A^, to the 
symmetric state \s). However, this coupling appears only for nonidentical 
atoms. 

From the master equation ()42|). we find that under the condition = 
\/r\r^ the equation of motion for the population of the antisymmetric state 
I a) is given by [11] 

(ri - r^)^ ^ .. , , 

Paa = J. , J. Pee + ^^e [psa " Pas) 

J- 1 + -L 2 

iPea-Pae) • (184) 

The equation ()184|1 shows that the non-decaying antisymmetric state \a) can 
be populated by spontaneous emission from the upper state |e) and also by 
the coherent interaction with the state \s). The first condition is satisfied 
only when Fi 7^ while the other condition is satisfied only when Ac 7^ 0. 
Thus, the transfer of population to the state \a) from the upper state |e) and 
the symmetric state |s) does not appear when the atoms are identical, but 
is possible for nonidentical atoms. 

We illustrate this effect in Fig. [T31 where we plot the steady-state pop- 
ulation of the maximally entangled state \a) as a function of A/, for two 
different types of nonidentical atoms. In the first case the atoms have the 
same damping rates (Fi = F2) but different transition frequencies (A 7^ 0), 
while in the second case the atoms have the same frequencies (A = 0) but 
different damping rates (Fi 7^ F2). It is seen from Fig. that in both cases 
the antisymmetric state can be populated even if is not directly driven from 
the ground state. The population is transferred to \a) through the coherent 
interaction Ac which leaves the other excited states completely unpopulated. 
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Figure 13: The steady-state population of the maximally entangled antisym- 
metric state \a) for Q = 10ri,ni2 = lOFi and r2 = Fi, A = Fi (solid line), 
F2 = 2Fi, A = (dashed line). 

This is shown in Fig. [T31 where we plot the steady-state populations pss and 
Pee of the states \s) and |e). It is apparent from Fig.[T31that at A^, = — ^^12 the 
states \s) and |e) are not populated. However, the population is not entirely 
trapped in the antisymmetric state |a), but rather in a linear superposition 
of the antisymmetric and ground states. This is illustrated in Fig. where 
we plot the steady-state population paa for the same parameters as in Fig. [T31 
but different fl. Clearly, for a small fl the steady-state population paa ~ |, 
and the amount of the population increases with increasing Q. The popula- 
tion paa attains the maximum value paa ~ 1 for a very strong driving field. 

This result shows that we can relatively easily prepare two nonidentical 
atoms in the maximally entangled antisymmetric state. The closeness of the 
prepared state to the ideal one is measured by the fidelity F. Here F is equal 
to the obtained maximum population in the state \a). For i7 ^ F the fidelity 
of the prepared state is maximal, equal to 1. As we have already mentioned, 
the system has the advantage that the maximally entangled state \a) does 
not decay, i.e. is a decoherence-free state. 

7.2.3 Atom-cavity-field interaction 

There have been several proposals to generate the antisymmetric state \a) in 
a system of two identical atoms interacting with a single-mode cavity field. 
In this case, collective effects arise from the interaction between the atoms 
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Figure 14: The steady-state populations of the upper state |e) (sohd hne) 
and the symmetric state |s) (dashed hne) for r2 = Ti,Q = 10ri,fii2 = lOFi 
and A = Fi. 

induced by a strong couphng of the atoms to the cavity mode. An excited 
atom emits a photon into the cavity mode that is almost immediately ab- 
sorbed by the second atom. Plenio et al. |155j have considered a system of 
two atoms trapped inside an optical cavity and separated by a distance much 
larger than the optical wavelength. This allows for the selective excitation of 
only one of the atoms. In this scheme the generation of the antisymmetric 
state relies on the concept of conditional dynamics due to continuous ob- 
servation of the cavity field. If only one atom is excited and no photon is 
detected outside the cavity, the atoms are prepared in a dark state |156j . 
which is equivalent to the antisymmetric state \a). 

In earlier studies, Phoenix and Barnett |157j . Kudrayvtsev and Knight |158j 
and Cirac and Zoller |159j have analyzed two-atom Jaynes-Cummings models 
for a violation of Bell's inequality, and have shown that the atoms moving 
across a single-mode cavity can be prepared in the antisymmetric state via 
the interaction with the cavity field. In this scheme, the preparation of the 
antisymmetric state takes place in two steps. In the first step, one atom 
initially prepared in its excited state |ei) is sent through a single- mode cav- 
ity being in the vacuum state |0)c. During the interaction with the cavity 
mode, the atomic population undergoes the vacuum Rabi oscillations, and 
the interaction time was varied by selecting different atomic velocities. If the 
velocity of the atom is such that the interaction time of the atom with the 
cavity mode is equal to quarter of the vacuum Rabi oscillations, the state 
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Figure 15: The steady-state population of the antisymmetric state \a) for 
T2 = ri,fii2 = lOFi, A = Ti and different il: O = Ti (sohd hne), Q = 5ri 
(dashed line), fl = 20ri (dashed-dotted line). 



of the combined system, the atom plus the cavity mode, is a superposition 
state 

\aic) = ^i\e,)\0)c'\gi)\l)c) . (185) 

Hence, the state of the total system, two atoms plus the cavity mode, after 
the first atom has crossed the cavity is 

l^i) = ;^(|ei)|0)e-M|l)c)M . (186) 

If we now send the second atom, being in its ground state, with the selected 
velocity such that during the interaction with the cavity mode the atom 
undergoes half of the vacuum Rabi oscillation, the final state of the system 
becomes 

I*i2c) = ^(|ei)|0)eb2)-|^i)|0)c|e2)) 

= ^(|ei)|^2)-M|e2))|0)e=|a)|0)e. (187) 

Thus, the final state of the system is a product state of the atomic antisym- 
metric state I a) and the vacuum state of the cavity mode. In this scheme 
the cavity mode is left in the vacuum state which prevents the antisymmet- 
ric state from any noise of the cavity. The scheme to entangle two atoms 
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in a cavity, proposed by Cirac and ZoUer |159j . has recently been realized 
experimentally by Hagly et al. [33|. 

Gerry |l(jUj has proposed a similar method based on a dispersive interac- 
tion of the atoms with a cavity mode prepared in a coherent state The 
atoms enter the cavity in superposition states 

W2) = ^(|e2)-^b)) . (188) 

After passage of the second atom, the final state of the system is 

I^i2c) = ^{(|^7i)|^72) + |ei)|e2))|-a) 

+ ^{\el)\g2)-\9l)\e2))\a)} • (189) 

Thus, if the cavity field is measured and found in the state the atoms 
are in the antisymmetric state. If the cavity field is found in the state |— a), 
the atoms are in the entangled state 

I^i2(-«)) = ^(|^?i)|^?2) + |ei)|e2)) . (190) 

The state p9()|l is called as a two photon entangled state. In section IHl we will 
discuss another method of preparing the system in the two-photon entangled 
state based on the interaction of two atoms with a squeezed vacuum field. 



7.3 Entanglement of two distant atoms 

In the previous subsection, we have discussed different excitation processes 
which can prepare two atoms in the antisymmetric state. The analysis in- 
volved single mode cavities, but ignored spontaneous emission from the atoms 
and the cavity damping. Here, we will extend this analysis to include spon- 
taneous emission from the atoms and the cavity damping |161j . We will show 
that two atoms separated by an arbitrary distance ri2 and interacting with 
a strongly damped cavity mode can behave as the Dicke model even if there 
is no assumed interaction between the atoms. 

Consider two identical atoms separated by a large distance such that 
~ and Q12 ~ 0. The interatomic axis is oriented perpendicular to 
the direction of the cavity mode (cavity axis) which is driven by an external 
coherent laser field of the Rabi frequency Q. The atoms are coupled to 
the cavity mode with coupling constant g, and damped at the rate F by 
spontaneous emission to modes other than the privileged cavity mode. For 
simplicity, we assume that the frequencies of the laser field ul and the cavity 
mode Uc are both equal to the atomic transition frequency ujq. The master 
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equation for the density operator pac of the system of two atoms plus cavity 
field has the form 



dpac 
dt 



a + a\'f)ac 
1 
2 



--in 

2 

f-.J- J-'aPac (-.J- c-^cPac ; 



5* + aS~^, , 



(191) 



where 



LaPa 
Lcpa 



(^PacS^S^ + S^S^ Pac — 25*^ PacSi 

1=1 

a^apac + PacCt^a - ^apaco) , 



(192) 



are operators representing the damping of the atoms by spontaneous emission 
and of the field by cavity decay, respectively; a and d'^ are the cavity-mode 
annihilation and creation operators, = Sf" + 5*^ are collective atomic 
operators, and Fc is the cavity damping rate. 

To explore the dynamics of the atoms, we assume the "bad-cavity" limit 
of Tc ^ g ^ T. This enables us to adiabatically eliminate the cavity mode 
and obtain a master equation for the reduced density operator of the atoms. 
We make the unitary transformation 



where 



Pt = D{-ri)pacDiv) 



(193) 



(194) 



is the displacement operator, and rj = fl/Tc- 

The master equation for the transformed operator reduces to 



dpT 
dt 



-igr] [S+ + S , pt\ - -ig 
-^TLaPT - ^TcLcPt ■ 



5*0^ + aS~^, Pt 



(195) 



We now introduce the photon number representation for the density operator 
Pt with respect to the cavity mode 



Pt = Pmn\m){n\ 

m,n=0 



(196) 



where pmn are the density matrix elements in the basis of the photon number 
states of the cavity mode. Since the cavity mode is strongly damped, we 
can neglect populations of the highly excited cavity modes and limit the 
expansion to m, n = 1. Under this approximation, the master equation p95j] 
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leads to the following set of coupled equations of motion for the density 
matrix elements 



Poo = LpQQ - ^Ag (S+pio - PmS'^ + TcPu , 

pio = LpiQ - ^ig (S^poo - Pii'S'") - ^r^pio , 

Pii = Lpii - ^ig (S^poi - pio5'+) - TcPii , (197) 

where Lpij = \igr] [5"+ + S' , pij] - \TLaPij. 

We note that the field-matrix elements pmn are still operators with respect 
to the atoms. Moreover 

Poo + Pii = Tr^;' (pt) = p (198) 

is the reduced density operator of the atoms. 

For the case of a strong cavity damping the most populated state of the 
cavity field is the ground state |0), and then we can assume that the coherence 
pio changes slowly in time, so that we can take pio = 0. Hence, we find that 



-L r 



(s-poo-piiS-) . (199) 
Substituting this solution to poo and pn, we get 



Poo = -f'Poo + TcPii - ^ [S^S Poo + PooS^S - 2S^piiS ) , 

pn = Lpu-r,pu + ^{2S-pooS^-S-S+pi,-pnS-S+) .(200) 

Adding these two equations together and neglecting population of the state 
|1), we obtain the master equation for the reduced density operator of the 
atoms as 



dp 1 . 



S+ + S- 



ItLp 



'S+S-p + pS+S- -2S-pS+) . (201) 



2r 



c 



The first term in Eq. ()20H) describes the interaction of the atoms with the 
driving field of an effective Rabi frequency grj. The second term represents 
the usual damping of the atoms by spontaneous emission, whereas the last 
term describes the damping of the collective system with an effective damping 
rate g'^/Tc- If we choose the parameters such that the collective damping is 
much larger than the spontaneous rates of the single atoms, the second term 
in Eq. (j20ip can be ignored, and then the master equation (j20ip describes 
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the time evolution of the collective two-atom system. Thus, two indepen- 
dent atoms located inside a strongly damped one-mode cavity behave as a 
single collective small sample model (Dicke model) with the damping rate 
(yf^/Fc. This model, however, requires that the atoms are strongly coupled 
to the cavity mode {g ^ T) and are located inside the cavity such that the 
interatomic axis is perpendicular to the direction of the cavity mode and the 
driving field. 

7.4 Preparation of a superposition of the antisymmet- 
ric and the ground states 

In the section 17.2.21 we have shown that two nonidentical two-level atoms 
can be prepared in an arbitrary superposition of the maximally entangled 
antisymmetric state \a) and the ground state \g) 

m=l\a) + ^/l-m9) . (202) 

However, the preparation of the superposition state requires that the atoms 
have different transition frequencies. Recently, Beige et al. |162j have pro- 
posed a scheme in which the superposition state |$) can be prepared in 
a system of two identical atoms placed at fixed positions inside an optical 
cavity. 

Here, we discuss an alternative scheme where the superposition state |$) 
can be generated in two identical atoms driven in free space by a coherent 
laser field. This can happen when the atoms are in nonequivalent posi- 
tions in the driving field, i.e. the atoms experience different intensities and 
phases of the driving field. For a comparison, we first consider a specific 
geometry for the driving field, namely that the field is propagated perpen- 
dicularly to the atomic axis (fc^ ■ ri2 = 0). We find from Eq. that in 
this case the collective states are populated with the population distribution 
Pee = Paa < Pss- The populatiou distribution changes dramatically when the 
driving field propagates in directions different from perpendicular to the in- 
teratomic axis ,85^ .90, . In this situation the populations strongly depend 
on the interatomic separation and the detuning A^. This can produce the in- 
teresting modification that the collective states can be selectively populated. 
We show this by solving numerically the system of 15 equations for the den- 
sity matrix elements. The populations are plotted against the detuning Al 
in Fig. El for the laser field propagating in the direction of the interatomic 
axis. We see from Fig. UHl that the collective states \s) and |e) are populated 
at Aj;^ = and A^ = Q12. The antisjTiimetric state is significantly populated 
at Al = —Q12, and at this detuning the populations of the states \s) and 
|e) are close to zero. Since paa < 1, the population is distributed between 
the antisymmetric and the ground states, and therefore at A^ = —Q12 the 
system is in a superposition of the maximally entangled state \a) and the 
ground state \g). 
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Figure 16: The steady-state populations of the collective atomic states of 
two identical atoms as a function of A^^ for the driving field propagating in 
the direction of the interatomic axis, Q = 2.5T, ri2/X = 0.08 and p, ± fi2: 
Pee (solid line), Paa (dashed line), (dashed-dotted line). 



Turchette et al. ^] have recently realized experimentally a superposition 
state of the ground state and a non-maximally entangled antisymmetric state 
in two trapped ions. In the experiment two trapped barium ions were side- 
band cooled to their motional ground states. Transitions between the states 
of the ions were induced by Raman pulses using co-propagating lasers. The 
ions were at positions that experience different Rabi frequencies Vti and VL2 
of the laser fields. By preparing the initial motional ground state with one 
ion excited 161)1(72)10), and applying the laser fields for a time t, the following 
entangled state |\E' {t)) was created 



+ 



sin (fit) |^)|1) + 
cos Vtt — 1 



9? 




cos fit - 1) 1 



|ei)|c/2) 



(203) 



where fi^ = fif + fi^. 

For fit = TT the entangled state ()203j) reduces to a non-maximally entan- 
gled antisymmetric state 



fi2 



|ei)|^2) - 



fi2 



9i)\e2) 



|0) 



(204) 
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Franke et al. |163j have proposed to use the non-maximally entangled 
state ()204j) to demonstrate the intrinsic difference between quantum and 
classical information transfers. The difference arises from the different ways 
in which the probabilities occur and is particularly clear in terms of entangled 
states. 

8 Detection of the entangled states 

In this section we describe two possible methods for detection of entangled 
states of two interacting atoms. One is the observation of angular intensity 
distribution of the fluorescence field emitted by the system of two interacting 
atoms. The other is based on quantum interference in which one observes 
interference pattern of the emitted field. Beige et al. jlH] have proposed a 
scheme, based on the quantum Zeno effect, to observe a decoherence-free 
state in a system of two three-level atoms located inside an optical cavity. 
The two schemes discussed here involve two two-level atoms in free space. 

8.1 Angular fluorescence distribution 

It is well known that the fluorescence field emitted from a two-atom system 
exhibits strong directional properties [TTl H?! IHH 1149] . This property can be 
used to detect an internal state of two interacting atoms. To show this, we 
consider the fluorescence intensity, defined in Eq. (jHH), that in terms of the 
density matrix elements of the collective atomic system can be written as 

I {R,t^ = u{R) {{pee + Pss) + cos {kri2 cos 9)] 

+ {pee + Paa) [1 " COS (fcria COS 9)] 

+ i{psa - Pas) sin {kri2 COS 9)} , (205) 

where 9 is the angle between the observation direction R and the vector fi2. 

The first term in Eq. (j2(J5j) arises from the fluorescence emitted on the 
|e) — s> |s) — i> \g) transitions, which involve the symmetric state. The second 
term arises from the |e) \a) — > \g) transitions through the antisymmet- 
ric state. These two terms describe two different channels of transitions for 
which the angular distribution is proportional to [1 ± cos (/cri2 cos 6')]. The 
last term in Eq. ()2()5j) originates from interference between these two radia- 
tion channels. It is evident from Eq. ()205|) that the angular distribution of 
the fluorescence field depends on the population of the entangled states \s) 
and I a). Moreover, independent of the interatomic separation ri2, the anti- 
symmetric state does not radiate in the direction perpendicular to the atomic 
for 9 = 7r/2 the factor [1 — cos {kri2 cos 6*)] vanishes. In contrast, the 
symmetric state radiates in all directions. Hence, the spatial distribution of 
the fluorescence field is not spherical unless pss = Paa and then the angular 
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distribution is spherically symmetric independent of the interatomic separa- 
tion. Therefore, an asymmetry in the angular distribution of the fluorescence 
field would be a compelling evidence that the entangled states \s) and \a) 
are not equally populated. If the fluorescence is detected in the direction 
perpendicular to the interatomic axis the observed intensity (if any) would 
correspond to the fluorescence field emitted from the symmetrical state |s) 
and/or the upper state |e). On the other hand, if there is no fluorescence 
detected in the direction perpendicular to the atomic axis, the population 
is entirely in a superposition of the antisymmetric state \a) and the ground 
state \g). 

Guo and Yang |164t 1165] have analyzed spontaneous decay from two 
atoms initially prepared in an entangled state. They have shown that the 
time evolution of the population inversion, which is proportional to the to- 
tal radiation intensity (f7T|). depends on the degree of entanglement of the 
initial state of the system. In Sections 14.1.11 and 14.1.21 we have shown that 
in the case of two non-identical atoms the time evolution of the total radia- 
tion intensity / (t) can exhibit quantum beats which result from the presence 
of correlations between the symmetric and antisymmetric states. In fact, 
quantum beats are present only if initially the system is in a non-maximally 
entangled state, and no quantum beats are predicted for maximally entangled 
as well as unentangled states. 

8.2 Interference pattern with a dark center 

An alternative way to detect entangled states of a two-atom system is to ob- 
serve an interference pattern of the fluorescence field emitted in the direction 
R, not necessary perpendicular to the interatomic axis. 

This scheme is particularly useful for detection of the symmetric or the 
antisymmetric state. To show this, we consider the visibility in terms of the 
density matrix elements of the collective atomic system as 

V = ^fi^fr . (206) 

Pss + Paa + 2p ee 

This simple formula shows that the sign of V depends on the population 
difference between the symmetric and antisymmetric states. For pss > paa the 
visibility V is positive, and then the interference pattern exhibits a maximum 
(bright center), whereas for p^s < Paa the visibility V is negative and then 
there is a minimum (dark center). The optimum positive (negative) value is 
V = 1 (V = — 1), and there is no interference pattern when V = 0. The later 
happens when p^s = Paa- 

Similar to the fluorescence intensity distribution, the visibility can provide 
an information about the entangled states of a two-atom system. When the 
system is prepared in the antisymmetric state or in a superposition of the 
antisymmetric and the ground states, pss = Pee = 0, and then the visibility 
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Figure 17: The visibility as a function of for ri2 = 0.1A,i7 = O.SF 
and various angles 6*^; 6^ = n/2 (solid line), 9l = t^/^ (dashed line), 9l = 
(dashed-dotted line). 

has the optimum negative value V = —1. On the other hand, when the 
system is prepared in the symmetric state or in a linear superposition of the 
symmetric and ground states, the visibility has the maximum positive value 
V = 1. 

The earliest theoretical studies of the fringe visibility involved a coherent 
driving field which produces an interference pattern with a bright center. Re- 
cently, Meyer and Yeoman jl66j have shown that in contrast to the coherent 
excitation, the incoherent field produces an interference pattern with a dark 
center. Dung and Ujihara |149j have shown that the fringe contrast factor 
can be negative for spontaneous emission from two undriven atoms, with ini- 
tially one atom excited. Interference pattern with a dark center can also be 
obtained with a coherent driving field |147j . This happens when the atoms 
experience different phases and/or intensities of the driving field. To show 
this, we solve numerically the master equation for the steady-state den- 
sity matrix elements of the driven system of two atoms. The visibility V is 
plotted in Fig. [T7|as a function of the detuning for = O.IA, Q = 0.25F 
and various angles 6l between the interatomic axis and the direction of prop- 
agation of the laser field. The visibility V is positive for most values of A^,, 
except A^ ^ —^12- At this detuning the parameter V is negative and reaches 
the optimum negative value V = —1 indicating that the system produces in- 
terference pattern with a dark center. In Fig. we plot the populations 
of the symmetric and antisymmetric states for the same parameters as in 
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Figure 18: Populations of the symmetric and antisymmetric states for the 
same parameters as in Fig. [T3 with Q]^ = 0. 

Fig 1171 It is evident from Fig. ITHl that at A^, = —Qu the antisymmetric 
state is significantly populated, whereas the population of the symmetric 
state is close to zero. This population difference leads to negative values of 
V, as predicted by Eq. ()206|1 and seen in Fig. El Experimental observation 
of the interference pattern with a dark center would be an interesting demon- 
stration of the controled excitation of a two-atom system to the entangled 
antisymmetric state. 

9 Two-photon entangled states 

In our discussions to date on entanglement creation in two-atom systems, we 
have focused on different methods of creating entangled states of the form 



As we have shown in Sec. 13.11 entangled states of the above form are generated 
by the dipole-dipole interaction between the atoms and the preparation of 
these states is sensitive to the difference A between the atomic transition 
frequencies and to the atomic decay rates. These states are better known as 
the symmetric and antisymmetric collective atomic states. 

Apart from the symmetric and antisymmetric states, there are two other 
collective states of the two-atom system: the ground state \g) = |5'i)|5'2) and 
the upper state |e) = |ei)|e2), which are also product states of the individual 




(207) 
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atomic states. These states are not affected by tlie dipole-dipole interaction 

^12- 

In this section, we discuss a method of creating entanglement between 
these two states of the general form 

\T) = Cg\g)±c,\e) , (208) 

where Cg and Cg are transformation parameters such that |cgp + |cep = 1. 
The entangled states of the form (|2U8|) are known in literature as pairwise 
atomic states [211 1221 12S1 12Z| or multi-atom squeezed states [23j. According 
to Eq. (jU^ . the collective ground and excited states are separated in energy 
by 2hujQ, and therefore we can call the states |T) as two-photon entangled 
(TPE) states. 

The two-photon entangled states cannot be generated by a simple co- 
herent excitation. A coherent field applied to the two-atom system couples 
to one-photon transitions. The problem is that coherent excitation popu- 
lates not only the upper state |e) but also the intermediate states \s) and 
I a), see Eq. fllSUp . The two-photon entangled states ()208|) are superpositions 
of the collective ground and excited states with no contribution from the 
intermediate collective states \s) and \a). 

The two-photon behavior of the entangled states (j2U8|) suggests that the 
simplest technique for generating the TPE states would be by applying a 
two-photon excitation process. An obvious candidate is a squeezed vacuum 
field which is characterized by strong two-photon correlations which would 
enable the transition \g) \e) to occur effectively in a single step without 
populating the intermediate states. We will illustrate this effect by ana- 
lyzing the populations of the collective atomic states of a two-atom system 
interacting with a squeezed vacuum field. 

9.1 Populations of the entangled states in a squeezed 
vacuum 

The general master equation (jH^ allows us to calculate the populations of 
the collective atomic states and coherences, which gives information about 
the stationary state of a two-atom system. We first consider a system of two 
identical atoms, separated by an arbitrary distance and interacting with a 
squeezed vacuum field. For simplicity, we assume that the carrier frequency 
Ug of the squeezed vacuum field is resonant to the atomic transition frequency 
uq, and the squeezed field is perfectly matched to the atoms, Diug) = 1 and 

es = 7i. 

From the master equation (jHBjl . we find the following equations of motion 
for the populations of the collective states and the two-photon coherences of 
the collective system of two identical atoms 

Pee = -2r (iV +l) Pee + iV [(F + Fia) Pss + (r - Fis) Paa] + ri2|M|p„ , 
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Pss = iT + Tu){N-{3N+l)pss-Np,a + Pee-\M\p^} , 
Paa = {T-Tu){N-{3N+l)paa~Npss+Pee+\M\p^} , 

Pu = 2ri2|M| - (2iv + i)rp„ 

-2\M\ [(r + 2Tu) Pss - (r - 2Tu) Paa] , (209) 

where N = N {uq), M = M{ujo) and pu = Peg exp(-i0j + pge exp{i(j)s). 

It is seen from Eq. fl209j) that the evolution of the populations depends 
on the two-photon coherencies peg and pge, which can transfer the population 
from the ground state \g) directly to the upper state |e) leaving the states 
|s) and I a) unpopulated. The evolution of the populations depends on ri2, 
but is completely independent of the dipole-dipole interaction Qi2- 

Similar to the interaction with the ordinary vacuum, discussed in Sec. l3.H 
the steady-state solution of Eqs. (j209p depends on whether = F or 7^ 
r. Assuming that ri2 = F and setting the left-hand side of equations ()209p 
equal to zero, we obtain the steady-state solutions for the populations and 
the two-photon coherence in the Dicke model. A straightforward algebraic 
manipulation of Eqs. (|209|) leads to the following steady-state solutions 



(2N + 1) 


1 - 1 


(2iV - 1) 


\M\ 


2 


(2N + 1) 


1 (sm + 3iv + 1 - 


3|M|2) 



iV(iV + l) - |M|2 

~ 3iV2 + 3iV + l-3|M|2 ' 

2|M| , , 

Pu = s , ~ — ^ ■ 210 

(2A^ + 1) (3A^2 + + 1 - 3|M|2J 

The steady-state populations depend strongly on the squeezing correla- 
tions M. For a classical squeezed field with the maximal correlations M = N 
the steady-state populations reduces to 

N 

Pss -- 



3N+1 

2iV2 

Pee = -1—^ \-r~' V • (211) 

{2N + 1) (3iV + 1) 

It is easily to check that in this case the populations obey a Boltzmann 
distribution with pgg > ps^ > pee- Moreover, in the limit of low intensities 
(A^ ^ 1) of the field, the population pee is proportional to A^^, showing 
that in classical fields the population exhibits a quadratic dependence on the 
intensity. 

The population distribution is qualitatively different for a quantum squeezed 
field with the maximal correlations |M|2 = A^(A^ + 1). In this case, the sta- 



76 



tionary populations of the excited collective states are 



N 



Pe 



(2N + l) 



(212) 



Clearly, the symmetric state is not populated. In this case the populations 
no longer obey the Boltzmann distribution. The population is distributed 
only between the ground state \g) and the upper state |e). Moreover, it can 
be seen from Eq. (|212|) that for a weak quantum squeezed field the popula- 
tion Pee depends linearly on the intensity. This distinctive feature reflects the 
direct modifications of the two-photon absorption that the nonclassical pho- 
ton correlations enable the transition \g) — * |e) to occurs in a "single step" 
proportional to N . In other words, the nonclassical two-photon correlations 
entangle the ground state \g) and the upper state |e) with no contribution 
from the symmetric state \s). 

The question we are interested in concerns the final state of the system 
and its purity. To answer this question, we apply Eq. ()210|) and find that in 
the steady-state, the density matrix of the system is given by 



P99 






Pss 





Pge 




(213) 



Peg U Pe, 

where pij are the non-zero steady-state density matrix elements. 

It is evident from Eq. (j213|) that in the squeezed vacuum the density 
matrix of the system is not diagonal due to the presence of the two-photon 
coherencies pge and p^g- This indicates that the collective states \g), \s) 
and |e) are no longer eigenstates of the system. The density matrix can be 
rediagonalized by including peg and pge to give the new (entangled) states 



[{Pi - pee) \g) + Peg\e)] / [(Pi - peef + \peg? 



[Pge\g) + {P2 

\s). 



Pgg)\<^)]l {P2- Pggf + \peg\' 



(214) 



where the diagonal probabilities are 



Pi 

P2 
P3 

In view of Eqs. 



1 

2 
1 

2 

Ps 



(Pgg + Pee) 



+ 



{pgg + Pee) - 



{pgg 
{pgg 



peef +^Peg[ 



peef + 4:\peg\^ 



(215) 



and 



it is easy to see that the squeezed vacuum 
causes the system to decay into entangled states which are linear superpo- 
sitions of the collective ground state \g) and the upper state |e). The in- 
termediate symmetric state remains unchanged under the squeezed vacuum 
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excitation. In general, the states ()214|) are mixed states. However, for perfect 
correlations |Mp = A^(A^+ 1) the populations P2 and P3 are zero leaving the 
population only in the state |Ti). Hence, the state |Ti) is a pure state of the 
system of two atoms driven by a squeezed vacuum field. From Eqs. (I214|) . 
we find that the pure entangled state |Ti) is given by |167j 



2iV + 1 



iV + l|^) + ViV|e) 



(216) 



The pure state ()216|) is non-maximally entangled state, it reduces to a max- 
imally entangled state for iV S> 1. The entangled state is analogous to the 
pairwise atomic state PH 123 1211 12Z| or the multi-atom squeezed state j2S] , 
(see also Ref. (221), predicted in the small sample model of two coupled atoms. 

The pure entangled state |Ti) is characteristic not only of the two- 
atom Dicke model, but in general of the Dicke model of an even number 
of atoms |168j . The A^-atom Dicke system interacting with a squeezed vac- 
uum can decay to a state which the density operator is given by 



p = Cn{lxS- + uS+) [fiS+ + uS-) , (217) 

if is odd, or 

p=|T)(T|, (218) 

if N is even, where Cn is the normalization constant, are the collective 
atomic operators, /i^ = z/^ + 1 = iV + 1, and |T) is defined by 

fiS- + uS^)\T) = . (219) 



Thus, for an even number of atoms the stationary state of the system is the 
pure pairwise atomic state. 



9.2 Effect of the antisymmetric state on the purity of 
the system 

The pure entangled state |Ti) can be obtained for perfect matching of the 
squeezed modes to the atoms and interatomic separations much smaller 
than the optical wavelength. To achieve perfect matching, it is necessary 
to squeeze of all the modes to which the atoms are coupled. That is, the 
squeezed modes must occupy the whole in solid angle of the space surround- 
ing the atoms. This is not possible to achieve with the present experiments 
in free space, and in order to avoid the difficulty cavity environments have 
been suggested EH] . Inside a cavity the atoms interact strongly only with 
the privileged cavity modes. By the squeezing of these cavity modes, which 
occupy only a small solid angle about the cavity axis, it would be possible 
to achieve perfect matching of the squeezed field to the atoms. 
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There is, however, the practical problem to fulfil the second requirement 
that interatomic separations should be much smaller than the resonant wave- 
length. This assumption may prove difficult in experimental realization as 
with the present techniques two atoms can be trapped within distances of 
the order of a resonant wavelength |^ Ej. As we have shown in 

Sec. 13.11 the dynamics of such a system involve the antisymmetric state and 
are significantly different from the dynamics of the Dicke model. 

For two atoms separated by an arbitrary distance ri2, the collective damp- 
ing 7^ r, and then the steady-state solutions of Eqs. (j2U9|) are 



iV2 a^\M\'^{AN + l 
2 + ^ 



G 



(2N +1 
NiN + l) a\M\A2{2N + l 



{2N + 1 



G 



Pa 



n(n + i) 

2N + 1 



2 (2N + lf + a 



G 



2a[2N +1]^ \M\ 



G 



where a = T12/T, and 



(220) 



G = (2iV+ 1)^ |(2iV+ l)V4|M|2 a^-(2N+l 



(221) 



This result shows that the antisymmetric state is populated in the steady- 
state even for small interatomic separations (ri2 ~ F). For large interatomic 
separations F12 ~ 0, and then the symmetric and antisjnnmetric states are 
equally populated. When the interatomic separation decreases, the pop- 
ulation of the state \a) increases, whereas the population of the state \s) 
decreases and pss = for very small interatomic separations. These features 
are illustrated in Fig. [THl^a), where we plot the steady-state populations as a 
function of the interatomic separation for the maximally correlated quantum 
squeezed field. We see that the collective states are unequally populated 
and in the case of small ri2, the state \a) is the most populated state of the 
system, whereas the state \s) is not populated. In Fig. ITW b). we show the 
populations for the equivalent maximally correlated classical squeezed field, 
and in this case all states are populated independent of ri2. 

This fact can lead to a destruction of the purity of the stationary state 
of the system. To show this, we calculate the quantity 

Tr (p^) = pI + pI + pI + pI + \pX , (222) 



79 



(a) 



(b) 




Figure 19: The steady-state populations of the collective atomic states as a 
function of for (a) quantum squeezed field with |Mp = A^(A^ + 1), (b) 
classical squeezed field with \M\ = N, and = 0.05, fl ± fi2: Pee (solid 
line), paa (dashed line), pss (dashed-dotted line). 

which determines the purity of the system. Tr(/3^) = 1 corresponds to a 
pure state of the system, while Tr(/3^) < 1 corresponds to a mixed state. 
Tr(p^) = 1/4 describes a completely mixed state of the system. In Fig. I20| 
we display Tr(p^) as a function of the interatomic separation ri2 for perfect 
correlations |Mp = N{N + 1) and various A^. Clearly, the system is in a 
mixed state independent of the interatomic separation. Moreover, the purity 
decreases as increases. 

For small interatomic separation, the mixed state of the system is com- 
posed of two states: the TPE state |Ti) and the antisymmetric state \a). We 
illustrate this by diagonalizing the steady-state density matrix of the system 



P 



f Pgg pge \ 

Paa 

p,, 

V Peg Pee / 



(223) 



from which we find the new (entangled) states 



= [(Pi -Pee) 1^) +Peg|e)]/ [(Pi -pee)V IPegl' 
IT2) = [Pge\g) + {P2 - Pgg) \e)] / \{P2 - PggY + \Peg\ 
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Figure 20: Tr(p^) as a function of the interatomic separation for |Mp = 
N{N + 1), /i ± fi2 and different N: N = 0.05 (solid line), N = 0.5 (dashed 
line), N = 5 (dashed-dotted line). 



1^) , 



(224) 



where the diagonal probabilities (populations of the entangled states) are 
Pi 



1 1 r 2 ^ 

- [Pgg + pee) + - [{Pgg " Pee) + 4 \peg\' 



1 1 r 2 

^^2 = - (Pgg + Pee) " - [{Pgg " Pee) + 4 \peg\ 

P3 = Pss , 
-P4 = Paa • 



(225) 



Note, that the states |Ti), IT2) and IT3) are the same as for the small 
sample model, discussed in the preceding section. This means that the pres- 
ence of the antisymmetric state does not affect the two-photon entangled 
states, but it can affect the population distribution between the states and 
the purity of the system. In Fig. we plot the populations Pi of the states 
|Tj) as a function of the interatomic separation. The figure demonstrates 
that in the case of a quantum squeezed field the atoms are driven into a 
mixed state composed of only two entangled states |Ti) and |a), and there 
is a vanishing probability that the system is in the states IT2) and In 
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Figure 21: Populations of the entangled states ()224|) as a function of the 
interatomic separation for (a) quantum squeezed field with \M\^ = N{N + 1), 
(b) classical squeezed field with \M\ = N, and /i ± fi2, N = 0.5. In both 
frames Pi (solid line), P2 (dashed line), P3 (dashed-dotted line), P4 (dotted 
line) . 

contrast, for a classical squeezed field, shown in Fig. l^lT b). the atoms are 
driven to a mixed state composed of all the entangled states. 

Following the discussion presented in Sec. 13. H we can argue that the 
system can decay to the pure TPE state |Ti) with the interatomic separation 
included. This can happen when the observation time is shorter than F"^. 
The antisymmetric state \a) decays on a time scale ~ (F — Fi2)~^, and for 
F12 ~ F the decay rate of the antisymmetric state is much longer than F~^. 
By contrast, the state |s) decays on a time scale ~ (F + Fi2)~^ , which for 
F12 ~ F is shorter than F~^. Clearly, for observation times shorter than 
F~^, the antisymmetric state does not participate in the interaction and the 
system reaches the steady-state only between the triplet states. Thus, for 
perfect matching of the squeezed modes to the atoms the symmetric state is 
not populated and then the system is in the pure TPE state |Ti). 
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9.3 Two-photon entangled states for two non-identical 
atoms 

We now extend the analysis of the population distribution in a squeezed 
vacuum to the case of two nonidentical atoms. For two nonidentical atoms 
with A 7^ and Fi = r2 = F, the master equation leads to the following 
equations of motion for the density matrix elements 



Pe 



-2F [N + Ij Pee + N\r {pss + Paa) + F12 (p, 
+F12IMI (peg 

(F + Fi2e^^*) [N - (3iV + 1 j Pss - Npaa + p 

—T\M\ l^pg^e"*'^^^''"^^^*^'^''' _|_ p^^g«[2{ws-^i)t+0s] 
-F12IMI (p^^e-*[2('^=~'^«)*+'^=l + pg^e'^2{^'~'^o)t+M 

, - Npss + Pe 



pa 



JAt 



g il2{u}s-u}o)t+(f)s] _|_ p ^il2{u}s-u}o)t+(f)s] 

'n - {3N + 1) p 



JAt 



N 



{3N + 1 



[2(t<Js-<^l)t + 0s] 



Peg 



F - Fi2e 
+F|M| (pege~* 
-F12IMI (p^^e-*[2(^=~^«)*+'^=l + p 

(Pge)* = Fi2|M|e*[2("=""^)*+'^=l - 
-F|M|e*PK-c.,)t+0.] (p^^-p^J 
-2Fi2|M|e*l'("^-"«)*+'^=] {pss + Pc 



J[2{LJs-l^0)t+<t>s 



2Ar+l Fp, 



(226) 



where tuo = | (cui + co'2). 

Equations ()226p contain time-dependent terms which oscillate at frequen- 
cies exp(±iAt) and exp[±2i(ti;s —u!o)t + (j)s]. If we tune the squeezed vacuum 
field to the middle of the frequency difference between the atomic frequen- 
cies, i.e. Us = {uJi +iJ2)/2, the terms proportional to exp[±2z(co's — LJo)t + (f)s] 
become stationary in time. None of the other time dependent components is 
resonant with the frequency of the squeezed vacuum field. Consequently, for 
A ^ F, the time-dependent components oscillate rapidly in time and average 
to zero over long times. Therefore, we can make a secular approximation in 
which we ignore the rapidly oscillating terms, and find the following steady- 
state solutions 



1 (2A^- 1 



4 2A^ + 1 




(2iV + l) -4a2|M|^ 



(2iV + l) -4a2|M|^ 
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Equations ()227p are quite different from Eqs. ()220|) and sliow tliat in the case 
of non-identical atoms tlie symmetric and antisymmetric states are equally 
populated independent of the interatomic separation. These are, however, 
some similarities to the steady-state solutions of the Dicke model that for 
small interatomic separations pss = Paa ~ 0, and then only the collective 
ground and the upper states are populated. 

To conclude this section, we point out that by employing two spatially 
separated non-identical atoms of significantly different transition frequencies 
(A ^ r), it is possible to achieve the pure TPE state with the antisymmetric 
state fully participating in the interaction. 

10 Mapping of entangled states of light on 
atoms 

The generation of the pure TPE state is an example of mapping of a state 
of quantum correlated light onto an atomic system. The two-photon correla- 
tions contained in the squeezed vacuum field can be completely transferred 
to the atomic system and can be measured, for example, by detecting fluc- 
tuations of the fluorescence field emitted by the atomic system. Squeezing 
in the fluorescence field is proportional to the squeezing in the atomic dipole 
operators (spin squeezing) which, on the other hand, can be found from the 
steady-state solutions for the density matrix elements. 

10.1 Mapping of photon correlations 

Equation ()227|) shows that the collective damping parameter plays the 
role of a degree of the correlation transfer from the squeezed vacuum to 
the atomic system. For large interatomic separations, ~ 0, and there 
is no transfer of the correlations to the system. In contrast, for very small 
separations, ~ F, and then the correlations are completely transferred to 
the atomic system. 

However, the complete transfer of the correlations does not necessary 
mean that the two-photon correlations are stored in the pure TPE state. 
This happens only for two nonidentical atoms in the Dicke model, where 
the steady-state is the pure TPE state. For identical atoms separated by 
a finite distance ru only a part of the correlations can be stored in the 
antisymmetric state. This can be shown, for example, by calculating of the 
interference pattern of the fluorescence field emitted by the system. Using 
the steady-state solutions ()220|) . we find that the visibility in the interference 
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Figure 22: The visibility V as a function of the interatomic separation ri2 for 
(a) a quantum squeezed field with |Mp = A^(A^ + 1), (b) a classical squeezed 
field with \M\ = N, fi± fi2 and different N: N = 0.05 (sohd line), N = 0.5 
(dashed line), = 5 (dashed-dotted hne). 



pattern is given by jl7()j 

2a|Mp 



V 



iV(2iV+ l)%2|M|2 



+ {2N + l) - (2iV + 1 



(228) 



The visibility is negative indicating that the squeezing correlations are 
mostly stored in the antisymmetric state. In Fig. |22l we plot the visibility V 
as a function of the interatomic separation for a quantum squeezed field with 
|M|2 = iV(iV + 1), Fig. E2ta), and a classical squeezed field with \M\ = N, 
Fig. |221^b). An interference pattern with a dark center is observed for small 
interatomic separations (ri2 < A/2) and with the quantum squeezed field the 
visibility attains the maximal negative value of V ~ —0.7 for ri2 < 0.3A. 
According to Eq. ()222p . at these interatomic separations the antisymmet- 
ric state is the most populated state of the system. The value V = —0.7 
compared to the possible negative value V = — 1 indicates that 70% of the 
squeezing correlations are stored in the antisymmetric state. In Fig. |22^b), 
we show the visibility for a classical squeezed field with \M\ = N. The visi- 
bility is much smaller than that in the quantum squeezed field and vanishes 
when N ^ oo. In contrast, for the quantum squeezed field V approaches 
— 1/2 when N ^ oo. Thus, the visibility can provide the information about 
the degree of nonclassical correlations stored in the entangled state \a). 
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10.2 Mapping of the field fiuctuations 

The fluctuations of tlie electric field are determined by the normally ordered 
variance of the field operators as 123 ISOl ISI] 

(: {AEef ■.) = Y.Eu {^a\a^us) + {hMs)^''' + («L«L)^"'^') " (229) 

ks 

Using the correlation functions ()17p of the three-dimensional squeezed vac- 
uum field and choosing 9 = 7t/2, the variance of the incident squeezed vacuum 
field can be written as 

: {AEl^/.f ■) = 2Eo{n~\M\) , (230) 



where Eq is a constant. Since \M\ = \JN{N + 1) > N, the variance ()230|) is 
negative indicating that the incident field is in a squeezed state. 

On the other hand, the normally ordered variance of the emitted fiuores- 
cence field can be expressed in terms of the density matrix elements of the 
two-atom system as 

(A^r*)' :^ = Eo (2p,, + 2pee + \Pu\ cos2^) . (231) 

Using the steady-state solutions ()210p and choosing 9 = tc/2, we find 

/ / ^2 \ (N- \M\) 

(: :) = 2E,^^^-^ . (232) 

Thus, at low intensities of the squeezed vacuum field (A^ ^ 1) the fiuctuations 
in the incident field are perfectly mapped onto the atomic system. For large 
intensities (A^ > 1), the thermal fiuctuations of the atomic dipoles dominate 
over the squeezed fiuctuations resulting in a reduction of squeezing in the 
fiuorescence field. 

The idea of mapping the field fiuctuations on the collective system of two 
atoms have been extended to multi-level atoms. For example, Kozhekin et 
al. |171j proposed a method of mapping of quantum states onto an atomic 
system based on the stimulated Raman absorption of propagating quantum 
light by a cloud of three-level atoms. Hald et al. jl72t I173j have experimen- 
tally observed the squeezed spin states of trapped three-level atoms in the 
A configuration, irradiated by a squeezed field. The observed squeezed spin 
states have been generated via entanglement exchange with the squeezed 
field that was completely absorbed by the atoms. The exchange process was, 
however, accomplished by spontaneous emission and only a limited amount 
of spin squeezing was achieved. Fleischhauer et al. |174j have considered a 
similar system of three-level atoms and have found that quantum states of 
single-photon fields can be mapped onto entangled states of the field and 
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the collective states of the atoms. This effect arises from a substantial reduc- 
tion of the group velocity of the field propagating through the atomic system, 
which results in a temporary storage of a quantum state of the field in atomic 
spins. These models are two examples of the continuing fruitful investigation 
of entanglement and reversible storage of information in collective atomic 
systems. 

11 Conclusions 

In this paper, we have reviewed the recent work on entanglement and non- 
classical effects in two-atom systems. We have discussed different schemes for 
generation of nonclassical states of light and preparation of two interacting 
atoms in specific entangled state. In particular, we have presented differ- 
ent methods of preparing two atoms in the antisymmetric state which is an 
example of a decoherence-free entangled state. The ability to prepare two- 
atom system in the decoherence-free state represents the ultimate quantum 
control of a physical system and opens the door for a number of applications 
ranging from quantum information, quantum computing to high-resolution 
spectroscopy. However, the practical implementation of entanglement in in- 
formation processing and quantum computation requires coherent manipu- 
lation of a large number of atoms, which is not an easy task. Although 
the two-atom systems, discussed in this review, are admittedly elementary 
models, they offers some advantages over the multiatom problem. Because 
of their simplicity, we have obtained detailed and almost exact solutions 
that can be easily interpreted physically, and thus provide insight into the 
behavior of more complicated multiatom systems. Moreover, many results 
discussed in this review is analogous to phenomena that one would expect 
in multiatom systems. For example, the nonexponential decay of the total 
radiation intensity from two nonidentical atoms is an elementary example of 
superradiant pulse formation, and a manifestation of the presence of coher- 
ences between the collective entangled states. A number of theoretical studies 
have been performed recently on entanglement and irreversible dynamics of 
a large number of atoms jll^ IT7T1 IT7^ ITTl ITTl IT7H1 ITTTl ITT^ ITTH^ 
These studies, however, have been limited to the Dicke model that ignores 
antisymmetric states of a multi-atom system. Nevertheless, the calculations 
have shown that population (information) can be stored in the collective 
atomic states or in the so called dark-state polaritons |18H I182j , which are 
quasiparticles associated with electromagnetically induced transparency in 
multi-level atomic systems. 
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